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ABSTRACT 

We investigate properties of stationary aligned and unaligned spiral perturbation con- 
figurations in a composite system of gravitationally coupled stellar and gaseous singu- 
lar isothermal disks (SIDs) using the two-fluid formalism. Both SIDs are taken to be 
razor thin and are in a self-consistent background equilibrium with power-law surface 
mass densities and flat rotation curves. We obtain stationary perturbation solutions 
for aligned and unaligned spiral logarithmic configurations in such a composite SID 
system and derive analytically existence criteria for these solutions. In comparison 
with the similar problem of a single SID studied by Shu et al. (2000), there are now 
two possible sets of solutions owing to an additional SID. For physically valid solutions, 
we explore parameter regimes involving the squared ratio [3 of velocity dispersions and 
the ratio i5 of the surface mass densities of the two SIDs. In terms of transition criteria 
from axisymmetric equilibria to aligned secular and spiral dynamical barlike instabil- 
ities, the corresponding T/|W| ratio of rotation to potential energies for a composite 
SID system depends on [3 and (5, varies in a wide range, and can be considerably 
lower than the oft-quoted value ~ 0.14. For both aligned and unaligned cases with 
azimuthal periodicities \m\ > 2, there exist certain parameter regimes where only one 
set of solutions is physically meaningful. For unaligned cases, we study marginal sta- 
bilities for axisymmetric (m = 0) and nonaxisymmetric (m ^ 0) disturbances. The 
resulting marginal instability curves, varying with parameters, are different from those 
of a single SID. The case of a composite partial SID system is also studied to include 
the gravitational influence of a dark-matter halo on the system equilibrium. For galac- 
tic applications, our model analysis contains more realistic elements and offers useful 
insights for the dynamics of disk galaxies consisting of stars and gas. Our analytical 
solutions are valuable for testing and benchmarking numerical codes. Starting from 
these solutions, numerical simulations are powerful to explore nonlinear dynamics such 
as large-scale spiral shocks. 

Key words: stars: formation — ISM: general — galaxies: kinematics and dynamics 
— galaxies: spiral — galaxies: structure. 



1 INTRODUCTION 



The theoretical model problem here is one that involves 
large-scale stationary density waves in a composite system 
of gravitationally coupled stellar and gaseous disks with the 
two fluid disks being both idealized as razor-thin singular 
isothermal disks (SIDs). For the gravitational effect of a 
background axisymmetric dark-matter halo, we also con- 
sider a background composite system of two coupled partial 
SIDs (Syer & Tremaine 1996; Shu, Laughlin, Lizano & Galli 



2000; Lou 2002). We search for stationary coplanar pertur- 
bation configurations in such a composite SID system. 

Over nearly four decades, there have been numerous 
theoretical studies on perturbation and stability properties 
of a composite system of stellar and gaseous disks, mostly in 
galactic contexts. The local dispersion relation for galactic 
spiral density waves in a composite disk system of stars and 
gas was first derived by Lin & Shu (1966, 1968), where the 
collective behavior of stars was described by a stellar dis- 
tribution function while the gaseous disk was treated as an 
isothermal fiuid. Their pioneer work was later followed up 
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by Kato (1972), who examined oscillations and overstability 
of density waves in a similar formalism. Using the two-fluid 
approach, Jog & Solomon (1984a, b) studied the growth of 
local a^xisymmetric perturbations in gravitationally coupled 
stellar and gaseous disks. They found that a composite disk 
system can be unstable owing to the gravitational coupling, 
even though the stellar and gaseous disks can be separately 
stable. Bertin & Romeo (1988) studied the role of gas on 
global spiral modes in a two-fluid model framework. Vander- 
voort (1991a,b) investigated effects of interstellar gas on os- 
cillations and stability of spheroidal galaxies. Romeo (1992) 
considered the stability of a two-component disk of finite 
thickness. The two-fluid formalism has been adopted into a 
modal analysis for morphologies of disk galaxies (Lowe et 
al. 1994). Elmegreen (1995) and Jog (1996) simultaneously 
approached a similar stability problem to derive an effective 
Qeff parameter (Safronov 1960; Toomre 1964) for axisym- 
metric two-fluid instabilities relevant to a disk galaxy. Lou 
& Fan (1998b) explored basic properties of open and tight- 
winding spiral density-wave modes in a composite disk sys- 
tem using the two-fluid formalism. 

Since Mestel (1963), the concept of SIDs has at- 
tracted considerable theoretical interests in various con- 
texts of disk dynamics in general (Zang 1976; Toomre 1977; 
Lemos, Kalnajs & Lynden-Bell 1991; Syer & Tremaine 1996; 
Lynden-Bell & Lemos 1999; Goodman & Evans 1999; Shu 
et al. 2000; Galli et al. 2001). These studies provide useful 
information for the research of star-formation (e.g. Shu et 
al. 1999 on formation and collapse of cloud cores in the birth 
of stars and planetary systems) and galactic structure com- 
munities (e.g., Bertin & Lin 1996 on the structure of barred 
and spiral galax:ies and Crane et al. 1993 on the light cusps 
seen in the nuclei of galaxies), and deepen our understand- 
ing for the dynamics of self-gravitating configurations that 
have a power-law density distribution. In particular, Shu et 
al. (2000) derived perturbation solutions and performed a 
stability analysis on an isopedically magnetized SID with 
a flat rotation curve. The background equilibrium surface 
mass density distribution was assumed to bear a power-law 
radial variation. They found both stationary aligned and 
unaligned logarithmic configurations* and offer physical in- 
terpretations for the marginal instability curves in a single 
SID. In contrast to their work on a single SID system, we 
are mainly interested in the situation of a disk system com- 
posed of a stellar SID and a gaseous SID. As it is more real- 
istic to consider large-scale dynamics of stellar and gaseous 
disks in a disk galaxy, the investigation on such a composite 
SID system can reveal more useful information. Motivated 
by such a prospect and combining our prior experience of 
treating composite disk system (Lou & Fan 1998b, 2000a, 
b), we search for both stationary aligned and unaligned or 

* Aligned disturbances involve distorted streamlines with the 
maximum and minimum radii at different radial locations lined 
up in the azimuth, while for unaligned or spiral disturbances, dis- 
torted streamlines with the maximum and minimum radii shifted 
in aaimuth at different radial locations (Kalnajs 1973). 



spiral logarithmic configurations in a composite SID sys- 
tem and discuss their stability properties (Safronov 1960; 

Toomre 1964; Goldreich & Lynden-BeU 1965; Ostriker & 
Peebles 1973; Elmegreen 1995; Jog 1996; Shu et al. 2000). 

Large-scale dynamics of the gas disk is a necessary com- 
ponent in the overall density-wave scenario. Moreover, most 
of the observational diagnostics involve processes of gaseous 
interstellar medium (ISM) on various sub-scales. In more 
realistic terms, numerical simulation experiments are indis- 
pensable for studying linear and nonlinear dynamical pro- 
cesses. For numerical code development in particular, the 
analytical solutions here not only offer important physical 
insights but also serve as valuable tools of benchmarking. 
Complementarily, numerical simulations starting from or 
based upon these stationary perturbation solutions can lead 
to insights for equilibrium, stability and nonlinear processes 
(e.g., spiral shocks to trigger star formation activities). Nu- 
merical simulations for clusters of galaxies and numerical 
simulations for galactic disk dynamics, although on totally 
different scales, are similar in several fundamental aspects, 
that is, both involve massive dark-matter halo, N-body grav- 
itational interaction (galaxies as mass points for a cluster 
versus stars as mass points for a disk galaxy), and gas dy- 
namics. The main difference lies in the geometry involved, 
namely, grossly spherical geometry for a typical galaxy clus- 
ter and a disk geometry for a spiral galax;y. Therefore, with 
proper adaptation, numerical codes designed for a galax:y 
cluster can be applied to a disk galaxy or vice versa. 

For spiral galax;ies, a more realistic model would involve 
a magnetized gas disk gravitationally coupled to a stellar 
disk with differential rotation in the presence of a massive 
dark-matter halo. For a single magnetized SID (MSID) with 
a coplanar magnetic field, we have recently derived solutions 
for stationary aligned and logarithmic spiral configurations 
(Lou 2002) and pointed out slow MSID configurations can 
persist in an extended radial range of a disk with flat rota- 
tion curve (Lou & Fan 2002) as in the case of NGC 6946 for 
a spiral pattern of interlaced optical and magnetic field arms 
(e.g.. Fan & Lou 1996; Lou & Fan 1998; Frick et al. 2000). 
The present hydrodynamic problem is not only interesting 
by itself, but also serves as a necessary step for construct- 
ing more realistic stationary configurations in a composite 
system of a stellar SID and an MSID. 

The paper is structured as follows. In Section 2, we 
present the basic two-fluid formalism for the composite sys- 
tem of SIDs coupled by self-gravity, derive equilibrium prop- 
erties for both SIDs, and obtain the linearized equations for 
small disturbances. Aligned and unaligned solutions for sta- 
tionary perturbations in a composite SID system and their 
stability properties arc studied and discussed in Section 3. 
We explore various parameter regimes for physical solutions 
and compare them with the results of the single SID case. 
The analysis is extended to a composite partial SID system 
in Section 4. We summarize computational results and anal- 
ysis in Section 5 and discuss potential galactic applications. 
Specific details are included in Appendices A— D. 
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2 TWO-FLUID SIDS FORMALISM 

For a model containing sufficient physics and for matliemat- 
ical simplicity, we adopt the two-fluid formalism for large- 
scale stationary aligned or unaligned disturbances in a back- 
ground rotational equilibrium with axisymmetry. In deal- 
ing with singularities and resonances, it would be physically 
more accurate to adopt the formalism of distribution func- 
tions (e.g. Lin & Shu 1966; Julian & Toomre 1966; Binney 
& Tremaine 1987). For the purpose of this study, such ir- 
regularities and resonances do not arise and the two-fluid 
approach will offer useful information for us to learn. In this 
section, we provide the basic equations for the two-fluid sys- 
tem, composed of a stellar disk and a gaseous disk. Given 
qualifications and assumptions, equilibrium properties of the 
stellar and gaseous SIDs with flat rotation curves (allowed 
to be different in a consistent manner) arc summarized. We 
derive coplanar perturbation equations in both stellar and 
gaseous SIDs, respectively. 

2.1 Two sets of coupled fluid equations 

For expediency, the two SIDs located at 2: = are treated 
as infinitcsimally thin, which arc sometimes referred to as 
razor-thin disks. For physical variables, we shall use super- 
scripts and/or subscripts s to denote the stellar disk and 
superscripts and/or subscripts g to denote the gaseous disk. 
The two razor-thin rotating disks in a composite system are 
modelled as two fluids coupled through the mutual gravita- 
tional interaction. In the present formulation of large-scale 
perturbations, we ignore nonideal effects such as viscosity, 
resistivity, and thermal diffusion etc. Then the two fully non- 
linear equation sets for the stellar and gaseous disks can be 
written out using cylindrical coordinates (r,<^, z) in the z = 
plane. For the stellar disk, we have 

^^Vi|:(rEV) + lA(ET) = 0, (1) 



dt r dr 



dt dr dip 



dt dr difi 



E^a^"^*"'^'^ d^ 



In parallel, we have for the gaseous disk 
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The coupling of the two sets of fluid equations is due to the 
gravitational potential through Poisson's integral 

-GT.{r',ip,t)r'dr' 



4>{ 



r,ip,t) = ^dtpj 



(7) 



[r'^ -|- f 2 _ 2rr' cos('i/' — 'p)]^^^ 
where E = E^ + E* is the total surface mass density. In 
equations (1) — (7), E** is the stellar surface mass density. 



u" is the radial component of the fluid velocity, j'' is the 
a— component of the specific angular momentum, and Os is 
the stellar velocity dispersion (or an effective "isothermal 
sound speed"), a^E'' stands for an effective pressure in the 
polytropic approximation, is the total gravitational po- 
tential expressed in terms of Poisson's integral. For physical 
variables of the gaseous disk, we simply replace the super- 
script s by g systematically. Here, we assume that the stellar 
and gaseous disks interact mainly through the mutual grav- 
itational coupling on large scales (Jog & Solomon 1984a,b; 
Bertin & Romeo 1988; Romeo 1992; Elmegreen 1995; Jog 
1996; Lou & Fan 1998b, 2000a,b). 

2.2 Properties of an axisymmetric equilibrium 

Before a coplanar perturbation analysis, one needs to adopt 
a background rotational equilibrium for the composite SID 
system consistent with equations (1) — (7). We assume ax- 
isymmetric background SIDs for both stellar and gaseous 
disks, with the same form of power-law surface mass densi- 
ties (E oc r~^) yet with different flat rotation curves. The 
equilibrium properties of the composite system can then be 
derived from the basic equations of Section 2.1. In applica- 
tions, the divergence towards r — > may be bypassed by 
introducing a gradual transition from disk to bulge or arti- 
ficial inner cut-outs. In theoretical analyses, such a r — » 
divergence poses a challenge of understanding the SID sta- 
bility properties (Zang 1976; Toomre 1977; Lynden-Bell & 
Lemos 1999; Evans & Read 1998; Goodman & Evans 1999; 
Shu et al. 2000). In particular, Shu et al. (2000) have made 
a systematic investigation on the problem in an attempt to 
summarize and clarify the relevant issues and to resolve the 
discrepancies between the disparate viewpoints. 

Using equations (2) and (5) for the background equilib- 
rium with Uq = Uq = 0, Q,s = jo/r'^, and fig = j^/r^, we 
readily obtain 

Eg = a?(l + D?)/[27rGr(l + 5)] , 

Eg = a^(l + I3^)5/[27rGr(l + 5)] , 

where 5 = Eg/Eg is the background surface mass density 
ratio. Furthermore, one can write 



(8) 



Os = UsDs/r , 

fig = ttgDg/r , 

Ks = {(2n,/r)[d(r^n«)/dr]}i''2 = y^Qs , 
Kg = {{2Qg/r)[d{r^ng)/dr]}'/^ = V2Qg , 



(9) 



(10) 



a?(D?-hl) =a^(£>^-M) , (11) 

where Os and Og are the mean angular rotation speeds of 

the stellar and gaseous SIDs, respectively, Ks and Kg are the 
corresponding epicyclic frequencies, and and Dg are the 
dimensionless parameters characterizing the level of rotation 
for the stellar and gaseous disks, respectively. The surface 
mass densities Es and Eg of both stellar and gaseous SIDs 
take the power-law form of oc r~^. Note that condition (11) 
(also valid for a composite system of two coupled partial 
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SIDs discussed later) is derived from the equilibrium radial 
momentum equations (2) and (5) using the polytropic ap- 
proximation for both SIDs. It implies an intimate relation 
among the two SID rotation speeds, the stellar velocity dis- 
persion and the gas sound speed. In this aspect, it is different 
from the usual proscription for a composite disk system (e.g., 
Jog & Solomon 1984a, b; Bertin & Romeo 1988; Elmegreen 
1995; Jog 1996; Lou & Fan 1998b). The two SID rotation 
speeds will be different as long as 7^ ag. This may induce 
streaming instabilities when the difference between the two 
SID rotation speeds is sufficiently large. 

Here, we introduce two useful parameters for a com- 
posite system of two gravitationally coupled SIDs. The first 
one is the SID surface mass density ratio S = E^/Eq. For 
late-type mature spiral galaxies, gas materials arc less than 
the stellar mass, that is, S < 1. However, for primordial 
disk galaxies, the gas materials exceed the stellar mass in 
general. Therefore, in our analysis and computations, both 
cases of (5 < 1 and S > 1 are considered. The second param- 
eter is /3 = ag/ttg for the square of the ratio of the velocity 
dispersion of stellar disk to the sound speed of gas disk. 
Typically, the stellar velocity dispersion is larger than the 
sound speed of gas disk. We therefore take > Og or /3 > 1. 
By "isothermal" , we mean constant and constant . For 
actual spiral galaxies, this represents a gross simplification. 
In the special situation of = Og , it follows from condition 
(11) that _Df = Dg. Thus, the two SIDs may be treated as a 
single SID, because the stellar and gas disks rotate with the 
same speed. We will show later that this l3 = 1 case gives 
some familiar results of a single SID which can be regarded 
as the limiting case of a two-SID system. In our analysis, we 
shall only consider the case^ oi (3 > 1. 



2.3 Coplanar perturbation equations in SIDs 

We now generally consider small nonaxisymmetric pertur- 
bations, marked along a physical variable with a subscript 
1, in both stellar and gaseous SIDs. For example, 

S^=Eg + E?, E«=Eg + E?, (12) 



u" = Uq + ul, = Uq + uf , (13) 



r=jo+ji f=fo+jf, (14) 



E = E" E'' = (Eg -h Eg) + (E? + E?) 



(15) 



Substituting expressions (12) — (15) into full equations 
(1) — (7) and linearizing about the axisymmetric background 
denoted by a subscript 0, we derive 

9Ef 1 a 9E! Eg djt 



^ Background equilibrium condition (11) guarantees that > 



as long as > when /3 > 1. 



^0 



(17) 



dj 



1 , , ^ djl a / aEf \ 

dt + "2n:"i + = - Eg + ) (^^^ 



d 



for the stellar SID, and 



dt r dr 







(19) 
(20) 
(21) 



dt +'^20^"! + "'' 9^ - dipy-'T.i^''", 

for the gaseous SID, with the total gravitational potential 
perturbation given by 



0i(r,</j,t) = j)d^ 



-G(E? + T.fydr 



(22) 



/() + r'^ — 2rr' cos{^p — 'y?)]"''''^ 

Assuming a Fourier periodic form of exp[i(c<;t — nvp)] for 
perturbation solutions in general (after taking the real part), 
we write for coplanar perturbations in the stellar disk as 

El — //"(r) cxp[i{cjt — rmp)] , 

u'l — (r) cxp[i{ijt — rriip)] , (23) 

ji = J\r) exp[i{ujt - mip)] , 
for coplanar perturbations in the gaseous disk as 

Ef = /Lt^(r) exp[i(a;t — rmp)] , 

uf = U^{r) cxp[i{u;t - rnip)] , (24) 

ji — (r) cxp[i{ujt — mip)] , 

and for the total gravitational potential perturbation as 

4>i = V{r) oxp[i{ujt — rmp)] (25) 

in the SID plane at ^; = 0. By substituting expressions (23) — 
(25) into equations (16) — (22), we derive for the stellar disk 



- mO,)/i'' + —{rEoU") - imEg^ = , (26) 



i{u> — m^ls)U" — 2^ls 



i{u,-mQs)r + r^U' =im{a^,^+V] , (28) 



for the gaseous disk 



i{oj - mQg)n^ + i^(rEg(7«) - imEg^ = , (29) 
r dr r-^ 



JO 



i(c^ - mng)C/« - 2ng— = -- ( + 1/ ) , (30) 



i(a;-mag)J«+r^i7«=imfa^|g +y) , (31) 



'9 \ ^0 

and for the total gravitational potential perturbation 

-G{iJL' + iiO)r'dr' 



V{r) 



id^ 

J Jo [r'2 + r 



(32) 



-)- 7.2 _ 2rr' cosV']^''^ 
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We now use equations (27) and (28) to express U" and J" in 
terms of = ailJ-"/'^o + V for the stellar SID and similarly, 
use equations (30) and (31) to express and in terms 
of = ag/i»/Eg + V for the gaseous SID. The resulting 
expressions then become 



U' = 
and 

- 

r 



(w — mf2s)2 — /tf 



1 



- 2Qs — + {ui — mQs)-T- 
r dr 



{oj — mfis)^ — 
for the stellar SID, and 

i 



{oj — mfis)- 



U' = 



{CO- 



mfig)2 



207 dr 



d 



-2n„ |-(aJ — mfig) — 

r dr 



and 

^ 1 

r (w — mf2g)2 — 



(w — mf2„)— • 
r 



Kg d 
2n^dr 



(33) 



(34) 



(35) 



^'^ (36) 



for the gaseous SID, respectively. 

Substitution of expressions (33) and (34) into equation 
(26) leads to 



= (a) — mQs)tJ-^ + 



r dr 



(a;. 



rSg 



2^3 h (w 

r 



dr 



r[{ijj — mfls)^ — fvi] 



(37) 



for the stellar SID. Similarly, substitution of expressions (35) 
and (36) into equation (29) leads to 



= (uj — mVln)ijF + - 

r dr 



x|r- 



mQ.g)'- 
mEg 



„_ m , \ d 

2\l„ \- (ui — m\l„) — 

r dr 



r[{LU - mQ.gy - K?g\ 



{u! — mO,g) 



m 



$9 
*9 



(38) 



for the gaseous SID. Equations (37) and (38) are to be solved 
with Poisson's integral (32). 

For stationary solutions [ui = 0) with zero pattern 
speed, equations (37) and (38) above can be cast into the 
forms of equations (39) and (40) below by invoking the back- 
ground equilibrium conditions (8) — (11). That is. 



D'i{m2 



_(rr^ 
2)\r 



dr 



dr^ 



rii + 



l+Dj V 
27rG 1 + 5 



(39) 



for the stellar SID, and 



L>2(m2 



dr 



dr^ 



„ 1 + VS 







(40) 



for the gaseous SID. Again, equations (39) and (40) are to be 
solved simultaneously together with Poisson's integral (32). 



3 ALIGNED AND UNALIGNED CASES 

3.1 Aligned perturbation configurations 

Let us now obtain the aligned stationary density wave pat- 
terns from equations (32), (39) and (40). We note in partic- 
ular that aligned perturbations relate to purely azimuthal 
propagations of density waves (see Section 3.2 of Lou 2002). 

3.1.1 The |?n| = Case: axisymmetric disturbances 

We first examine the m = case. With u! = m = Q, the solu- 
tion to equations (26) — (32) takes the forms of U" = = 0, 
= Kt/r, At" = Kf/r, J" = K^r, J" = Kir, V = Klnr, 
where the ratios of constants Kf/K, K(/K, K2/K, and 
K2/K arc chosen such that equations (27), (30), (32) can 
be satisfied. However, such a "solution" merely represents a 
rescaling of one axisymmetric equilibrium to a neighboring 
axisymmetric equilibrium (Shu ot al. 2000). This rescaling 
is allowed by equations (32), (39) and (40) but uninteresting 
in the present context. We now turn to cases with \m\ > 1. 

3.1.2 Cases with \m\ > 1: nonaxisymmetric disturbances 

In power-law disks, we consider aligned perturbations that 
carry the same power-law dependence as the equilibrium 
SID does. By this assumption, we mean the following exact 
potential-density relation 



^Jis = 0-3 /r, 
Ms = f^g/r, 
V = -2TTGrHs/\m\ - 2nGrHg/\m\ 



(41) 



where cTg and (Tg are constant coefficients. It can be verified 
that Poisson's integral (32) is satisfied (Shu ot al. 2000). A 
direct substitution of equation (41) into equations (39) and 
(40) then leads to the following equations: 

'm^ „ d d^ 

H~ = ■ 



- 2^ - rj-^ ) (if^r// + Gar/i') , (43) 



where coefficients Hi, H2, Gi and G2 are defined by 



Hi = 



H2 = 



L>2(m2-2) 



L>|(m2 -2) 



1 - 



Dj + l 1 

|m| 1 + 5 

Dg + 1 5 

\m\ 1 + S 



Gi = 

G2 = 



Di + i 



£)i(m2-2)|m| 1 + 5 ' 

5 



Dg + ^ 



(44) 



(45) 



D2(m2 -2)|m| 1 + 5 



(1 — Him = Gim /x' 



2, ,9 



By substituting the forms of and given by equation 
(41) into equations (42) and (43) above, we readily obtain 

(46) 
(47) 



and 



(1 — H2m^)ijF = Gam^/i" 
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It then follows from equations (46) and (47) that 

{1- Him'^){l- H2m'^) = GiG2m'^ . (48) 

In the above procedure, we have assumed 5 ^ in order 
to eliminate /i" and /i* from both sides of equation (48). 
Otherwise the problem would reduce to that of a single SID. 

Using definitions (44) and (45) for the expressions of 
Hi, H2, Gi and G2, we obtain explicitly 



Dlim 
1 - 



2) 



1 - 



D'i + 1 



Z)|(m2 -2) 



i\ l + S 



1 + 5 



(49) 



L>i(m2 -2) D|(m2 -2) (1 
which can be further simplified to 

Dl + l 1 



5y 



1 - 



rflf 2 



2) 



|m| 1 ^ 



{Dl + \){Dl + l) 



1 + 5 
ra^5 



(50) 



(1+5)2 • 

By condition (11) for the background equilibrium, we have 

D'I = Q{D1 + 1) - 1 with d = al/a^. Equation (50) can be 
cast into a quadratic equation in terms oi y = D^, namely 



where 



C2t/^ + Cij/ + Co = 



C2 = (3{rn - 2){rn + \m\ - 2) 



Ci = ( - 2m^ + 4 - 



2m 4 2 
r — 2m + 6m — 

1 + 



Co 



m 2 4 

+ 2m — m 

1 + 



2m5 Am 
1+5 ~ \+5 
2m^ 
1 + 5 

2m 



+ 



2m^ 
1 + 5 

-4 , 
m-'5 







1+5 l+S 



2m 2 '2m'' 4 
+ :; r — 2m — r + 2m 



(51) 
(51o) 

(516) 

1 + 5 - 1 + 5 • 
The physical meaning is that for aligned nonaxisymmetric 
stationary configurations to exist in a composite SID sys- 
tem, condition (50), or equivalently, condition (51) must be 
satisfied for proper values of given specified parameters 
m, 5 and (3. In essence, a purely azimuthal propagation of 
density wave in opposite direction relative to the SID rota- 
tion needs to be counterbalanced by disk rotation in order 
to appear stationary in an inertial frame of reference. There 
are two possible sets of density waves in a composite disk 
system (e.g., Lou & Fan 1998b). Therefore, there are two 
possible disk rotation speeds that may produce stationary 
perturbation configurations in general. 

Mathematically, one can show that equation (51) always 
has two real solutions for y = D^, except for the special 
case of m = 1. The two real solutions are not necessarily 
both physically valid due to the nonnegative requirement 
of > 0. Nevertheless, it can be proven later that there 
always exists at least one physical solution with y = > 0. 



For the aligned case of m = 1, it turns out that 
C2 = Ci = Co = by definitions (51o) — (51c) and 
— S„/E5 according to equation (46). Equation (51) 
can therefore be satisfied for arbitrary values of D^. This 
is quite similar to the m = 1 case of a single SID stud- 
ied by Shu ct al. (2000). While noting earlier concerns of 
spurious artifacts of the analytical method about this re- 
sult (Zang 1976; Toomre 1977), Shu et al (2000) expressed 
a strong belief that such aligned eccentric displacements are 
possible alternative states of equilibria for extended SIDs. 
Whether such m = 1 case represents a trivial translation 
of the origin of coordinates, we note for example that con- 
tours of surface mass density (background plus perturba- 
tion ones) are given by E = (Eor)/r + a cos tp/r =constant 
according to equation (41), while contours of surface mass 
density resulting from a small shift a along the a;— axis of 
the origin from a; = 0, ?/ = 0toa; = a, j/ = are given by 
E — [(Eor)/r](l + a cos (ys/r) =constant. It is then clear that 
m = 1 aligned perturbations are not trivial translation of the 
origin of coordinates for the present model problem.* On this 
ground, it seems plausible that stationary aligned eccentric 
|m| = 1 configurations are possible alternative equilibria 
for a composite SID system, at least mathematically. We 
note, however, for a single partial SID (Lou 2002) and for 
a composite system of two coupled partial SIDs (discussed 
later), such stationary aligned eccentric \m\ = 1 configu- 
rations are not allowed. In other words, stationary aligned 
eccentric |m| = 1 configurations might be possible for proto- 
stellar disks where daxk-matter halos axe not involved, but 
may not occur for disk galaxies where massive dark-matter 
halos are known to exist. 

We now study cases of m > 2. As noted earlier, in a 
real spiral galaxy, the stellar velocity dispersion Us is usually 
larger than the sound speed Ug of the gas disk^, that is, 
/3 > 1 in our model. Theoretically, the case of /? = 1 means 
that the stellar disk and gaseous disk have the same velocity 
dispersion, and thus the same rotation parameter Dg = Dg 
by condition (11). In some sense, we may then treat the two 
SIDs as one single SID with E = E" + E^ . It is expected that 
the solution should have something in common with a single 
SID system (Shu et al. 2000). We consider below the case 
ol (3 = 1 for analytical solutions and for insights of general 
solution properties. 

With /3 = 1, equations (50) or (51) can be reduced to 

(|m| - l)[L'?(m' - 2) - m'][(|m| + 2)dI - \m\] = (52) 
for arbitrary 5 values. Note that expressions (51o), (516) and 



t In analyzing stability properties of Maclaurin disks (Takahara 
1976; Smith 1979), the marginal stationary case of m = 1 rep- 
resents a trivial translation of the origin of coordinates involving 
a uniform velocity perturbation as shown in footnote 14 on p319 

of Binncy & Trcmainc (1987). Using tlic same analytical crite- 
rion here, one can demonstrate tliat contours of tiie surface mass 
density (see cqns 5-99 and 5-117 of Binncy & Trcmainc 1987) are 
equivalent to a small shift of the origin of coordinates. 
§ For example, one may talce Og = 7km and Os = 30km s~^ 
for a typical late-type spiral galaxy. 
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(51c) are independent of 5 when /3 = 1. We therefore have 
two solutions = m^/{m^ — 2) and = |m|/(|m| + 2) for 
\m\ > 2, and arbitrary values of for |mj — 1. Note that 
subsonic rotation Dj = |m|/(|m| + 2) for \rn\ > 2 is simply 
the result of a single SID (equation (27) of Sim et al. 2000), 
and the supersonic rotation Dj — ni" jiju' — 2) for |m| > 2 
is a novel feature of a composite SID system. 

For the phase relationship between the two surface mass 
perturbations /x" and /n®, we note that equation (46) can be 
written in the form of 

_ [D?(m2 -2) -m2](l + r5) 



Gim2 



\Tn\{DI 



(53) 



For supersonic rotations with Dj = m^/(m^ — 2), one simply 
has jJp / n" = —1, while for subsonic rotations with = 
\m\/{\m\ + 2), one has ^V/"^ = ^ (Lou & Fan 1998b). 

We have noted earlier (Lou & Fan 2002; Lou 2002) that 
stationary aligned nonaxisymmetric configurations in an in- 
ertial reference frame correspond to purely azimutlial prop- 
agations of density waves counterbalanced by SID rotation. 
Physically, /x® = — simply means that surface mass den- 
sity perturbations of stellar and gaseous SIDs arc completely 
out of phase to reduce the effect of gravity. This in turn im- 
plies a faster azimuthal density wave speed and thus requires 
a faster disk rotation (larger D1) to maintain a stationary 
configuration. The case of /j^ = 5jf means that mass distur- 
bance in gaseous disk is in phase with the mass disturbance 
in stellar disk, and their ratio is the same as the ratio of 
the background surface mass densities of the two SIDs. As 
the effect of self-gravity is enhanced, the azimuthal density 
waves speed is slower and thus a slower disk rotation (smaller 
DI) is needed to sustain a stationary configuration. 

For more realistic situation of /3 > 1 in general, we 
readily obtain two branches of solution for y = from 
quadratic equation (51) as functions of f5 when |m| > 2 and 5 
are specified. It is found that both solution branches (upper 
branch yi and lower branch y2) monotonically decrease for 
increasing values of /3, and asymptotically approach different 
limits when /3 goes to infinity (see three numerical examples 
shown in Fig. 1 and Appendix B for the variation trends). 
It can be shown analytically that the larger solution branch 
j/i (i.e., the higher branch) satisfies the following inequality 



2m(m -I- 1) 



< J/i < 



(54) 



m?-2 {rrfi -2){m + 2){l+5) ^ ^ m? -2 
where the lower bound on the left-hand side, which is always 
positive, is obtained by taking the limit of /3 ^ cx) and 
the upper bound on the right-hand side is determined by 
simply taking (3 = 1. One can also show analytically that the 
smaller solution branch j/2 (i-e., the lower branch) satisfies 
another inequality 

-1< 2/2 < \Tn\K\M +2) , (55) 
where the lower bound on the left-hand side is derived by 
letting /3 — » oo and the upper bound on the right is obtained 
by simply taking /3 = 1. Note that the left-hand side of yi 
in inequality (54) is always greater than the right-hand side 
of 2/2 in inequality (55), so that the two solution branches 
will not intersect with each other and yi remains always 
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Figure 1. Two sets of solution curves j/i and j/2 versus /3 for 
|m| = 2 and S = 0.2, 1, 5, respectively, where /3 varies in the 
interval [1, oo). The higher yi branch is always positive, while the 
lower j/2 branch becomes negative when /3 exceeds some critical 
value /3c defined by equation (56) (see also Table 1). Both solution 
branches decrease with increasing /? and approach different limits 
as /? — > oo. Note also that the yi branch approaches the limiting 
value much faster than the j/2 branch does. 



greater than j/2- In other words, for a given \m\, while there 
are two possible aligned solutions at maximum, there is only 
one aligned configuration that a composite SID system can 
support at a time. However, this docs not exclude the possi- 
bility that for different values of |mj, more than one aligned 
configurations may be sustained in a composite SID system 
simultaneously. 

By inequalities (54) and (55) and by numerical exam- 
ples shown in Figure 1, the lower branches are subsonic while 
the upper branches are positive but may vary from super- 
sonic (3 ^ 1) to subsonic for sufficiently large 13 values when 
6 is small enough (sec the case of m = 2 and 6 = 0.2). For 
late- type spiral galaxies, 6 ranges from 0.05 to 0.1 and (3 
may take values between ~ 15 — 20. Based on the variation 
trend displayed in Fig. 1, it is then possible for late- type 
disk galaxies to support stationary bar configurations of the 
upper branch with subsonic SID rotations. It is also inter- 
esting to infer that for early-type young disk galaxies with 
large values of 6, stationary bar configurations of the upper 
branch may be sustained by supersonic SID rotations. 

By examples of Fig. 1, it also becomes clear that for 
cases of |m| > 2, the lower solution branch 3/2 of equation 
(51) may become negative when l3 is larger than a critical 
value Pc for given values of m and 5. In other words, there is 
only one solution branch for (3 > j3c, i.e. the upper 3/1 = 
branch, which is physically meaningful; being negative, the 
lower j/2 = branch becomes unphysical and should be 
discarded. This critical value f3c can be determined analyti- 
cally in terms of 5 and m as 

_2{rn + 1) 



(m-l-2) 



1 + ' 



(56) 



(m? + 2m)5 + m? - 2 

For a given value of m, the critical value of f3c decreases 
with increasing values of 5. When 5 approaches infinity, the 
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critical value l3c goes to a limiting value 

/3cl™ = 2(m + l)/(m + 2) . (57) 

In contexts of disk galaxies, this seems to suggest that 
aligned barred configurations with m — 2 of the lower so- 
lution branch may be sustained by subsonic SID rotations 
only when /3 < 3/2 which is rather restrictive. 

Wc now examine the phase relationship between the 
surface mass density perturbations jj,^ and /i" in general. 
One can demonstrate that for the upper j/i = branch, 
the inequality 



/x« (rn^ + \rii\ - 2)6 + in'^ -2 ^ ' 

holds. It is found that jJ,^ / n" ratio decreases with increasing 
values of (see Appendix C) and therefore increases with 
increasing values of f3. The lower bound on the left-hand 
side of inequality (58) is determined by taking l3 = 1 and 
the upper bound on the right-hand side of inequality (58) is 
obtained by letting /3 — »• oo. For the lower j/2 ~ Di branch, 
we have the following inequality 

5 < Itx' <\m\{l + 5)-l (59) 

for perturbation surface mass density ratio As the 

mass ratio / ijf also increases with increasing values of (3, 
the lower bound on the left-hand side of inequality (59) is de- 
termined by taking 13 = 1. However, the upper bound on the 
right-hand side of inequality (59) is determined according to 
the possible maximum value of (3 that makes physically 
meaningful, i.e. the critical value pc given by expression (56). 
It is now clear that the = yi solution branch is always 
characterized by a negative fJ,^ /fJ-^ and the = j/2 solution 
branch is always characterized by a positive fi^ / . In this 
regard, the (3 = 1 case which can be analyzed thoroughly 
shows a general solution property. 

As observational diagnostics, our analysis suggests two 
different types of stationary aligned bar configurations in 
terms of surface mass density distributions. For stationary 
bar configurations of the upper branch solution, gas baxs 
or young stellar bars should be out of phase relative to old 
stellar bars. For stationary bar configurations of the lower 
branch solution, gas bars or young stellar bars should overlap 
with old stellar bars. Another point of interest is that for 
late-type spiral galaxies, f3 may be as large as ~ 15 — 20. 
By Fig.l, the lower y2 solution branches arc excluded and 
the secular barlike instabilities more likely occur along the 
upper j/2 solution branches. 

In general, the two ratios (3 and S are independent of 
each other. When 5 and (3 are specified, equation (51) can 
be solved for two values of Dg with different values of \m\. 
We now take on a specific numerical example below (see Fig. 
1). For \m\ = 2 and 5 = 5, equation (51) becomes 

Sfy'' + {-U-6l3)y- 14/3 + 22 = , (60) 

which has two solutions for = y, namely, 

yi = [7 3/3 -h (49 - 134/3 + 121/3^)i/^]/(8/3) (61) 
for -1 < n^/n" < -5/11, and 

y2 = [7 -h 3/3 - (49 - 134/3 + 121/3')'/']/(8/3) (62) 



for 5 < / < 11. For any value /3 > 1, j/i given by equa- 
tion (61) is always positive. For 1 < (3 < 11/7, j/2 given 
by equation (62) is nonncgativc, while when (3 > 11/7, y2 
becomes negative and thus unphysical. Besides other con- 
straints, this may imply that in a disk galaxy where the 
stellar velocity dispersion is much greater than the sound 
speed of the gaseous disk, the / < branch may tend 
to manifest. Note that the range of //®/a«^ is also fully de- 
termined by inequality (58) once |mj and 5 are known. 

For an numerical example of a late-type disk gala:xy, we 
follow the similar procedure shown above and take param- 
eters rn = 2, 5 = 0.05 and /3 = 4. Now equation (51) yields 
two solutions, namely, j/i = 0.6008 and j/2 = —0.2972. Ap- 
parently, the second j/2 solution should be discarded since 
P exceeds the critical value j3c given by equation (56), i.e., 
/3c = 11/4. The 2/1 solution is physically meaningful with 
Ijfi / = —0.0822. In this case, stationary aligned perturba- 
tions in gas surface mass density and stellar surface mass 
density are out of phase with each other. 

To identify the nature of the aligned solution condition 
(50) or (51), we now examine the closely relevant case stud- 
ied by Shu et al. (2000) and reemphasize the perspective 
that stationary aligned perturbation configurations should 
be regarded as purely azimuthal propagation of hydrody- 
namic density waves (Lou 2002; Lou & Fan 2002). For this 
purpose, we write solution condition (50) in a physically 
suggestive form of 

[n2(m^ - 2) - rnailr^ + 27rGS°|m|/r] 
X \^l(rn - 2) - -rnajr- + 27rGE°|m|/r] 
= 47r2G^E^EgmVr'' . (63) 

The right-hand side of equation (63) represents the mutual 
gravitational coupling between the two SIDs. In the absence 
of this coupling, the two factors on the left-hand side would 
give rise two separate conditions for stationary aligned per- 
turbation configurations with |m| > 2 for stellar and gaseous 
SIDs, respectively. For a single SID, be it stellar or gaseous, 
equation (63) therefore reduces to the form of 

m^n^ = K^-hm^o^/r^ -27rGEo|m|/r . (64) 

According to the well-known dispersion relation of density 

waves derived under the tight-winding or WKBJ approxi- 
mation (Lin & Shu 1964, 1966 or equation (39) of Shu et al. 
2000), equation (64) can be recovered by replacing the radial 
wavenumber |fc| with the azimuthal wavenumber |m|/r and 
setting w = in an inertia frame of reference as noted by Lou 
(2002) in the study of stationary MHD perturbation config- 
urations in a single MSID. By this procedure, it is quite 
clear that equation (64) describes an azimuthal propagation 
of hydrodynamic density waves, and a stationary pattern in 
the sidereal frame of reference requires specific values of 
for different |m| values (Lou 2002; Lou & Fan 2002). In ref- 
erence to dispersion relation (3.20) of Lou & Fan (1998b) 
with a proper adjustment of two different rotation speeds of 
SIDs, it is also transparent that equation (63) represents a 
purely azimuthal propagation of density waves in a compos- 
ite SID system. The two coupled SIDs rotate in such a way 
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|m| 


<5 


Critical value /3c for 
the lower 2/2 branch 


1 


2 


0.2 
1 

5 


2.3333 
1.8000 
1.5714 




oo 


1.5000 


3 


0.2 
1 

5 


2.0800 
1.8182 
1.6585 




oo 


1.6000 



Table 1. Critical value /3c of the lower j/2 branch solution for 
stationary aligned configurations in a composite SID system, de- 
termined by equation (56). See Fig. 1 for three examples. 

to render the azimuthal density wave pattern stationary in 
a sidereal frame of reference. 

In the context of a single SID, Shu et al. (2000) sug- 
gested that stationary condition (64) represents the onset 
of bifurcations of axisymmetric SIDs to nonaxisymmetri- 
cal SIDs that are more centrally condensed when the ro- 
tation rate D is systematically increased. This is analogous 
to classical Maclaurin-Jacobi or Maclaurin-Dedekind bifur- 
cations from spheroids to ellipsoids (Chandrasekhar 1969; 
Tassoul 1978) and involve secular instabilities that may be 
induced by viscous dissipations or by gravitational radiar 
tion (Bardeen et al. 1977). Shu et al. (2000) indicated that 
the instability mechanisms for aligned and spiral instabilities 
are fundamentally different and emphasized that wave prop- 
agation plays no role in aligned perturbations in contrast 
to spiral perturbations (see discussions in their subsection 
3.1). Base on our analysis, we would like to clarify here that 
both types of perturbations involve propagations of density 
waves, even though a radial wave propagation for the spiral 
case may lead to different instability mechanisms than those 
without radial wave propagation for the aligned case. It is 
also important to note that for stationary configurations to 
appear in a sidereal reference frame, these density waves 
travel in opposite sense of SID rotation relative to the SID 
system and may strike a balance between wave propagation 
and SID rotation. 

3.1.3 Secular Barlike Instabilities 

Strictly speaking, we have constructed analytically station- 
ary aligned nonaxisymmetric configurations in a composite 
SID system. Are they stable or do they merely represent 
transition states between axisymmetric equilibria and non- 
axisymmetric configurations? This is a challenging question. 
In the single SID case, it was suggested (Shu et al. 2000; Galli 
et al. 2001) that these solutions signal onsets of bifurcations 
from an axisymmetric SID to nonaxisymmetric SIDs (such 
as eccentric, oval, triangular distortions corresponding to 
m = 1, 2, 3 etc.). Based on transmisions and overreflections 



of leading and/or trailing spiral density waves across the 
corotation in a time-dependent problem, Shu et al. (2000) 
also made a novel suggestion that the condition for station- 
ary unaligned or logarithmic spiral configurations in a SID 
would determine whether spiral density waves may be swing- 
amplified (Goldrcich & Lyndon-Bell 1965; Fan & Lou 1997). 
In addition, Shu et al. (2000) examined the parameters of 
these stationarity conditions in reference to the disk stability 
criterion hypothesized by Ostriker & Peebles (1973) for the 
onset of bar- type instabilities (Miller et al. 1970; Hohl 1971; 
Kalnajs 1972) and estimated Ostriker-Peebles criterion for 
aligned secular and spiral dynamic barlike instabilities in a 
SID. Judging the similarity and difference between a single 
SID and a composite SID system, we examine stability prop- 
erties for stationary aligned and unaligned configurations in 
the same spirit. 

It should be noted that onsets of or transitions to bar- 
like instabilities may also occur at non-zero pattern speeds 
such as those of Maclaurin spheroids to Jacobi ellipsoids 
(Chandrasekhar 1969) or those of Kalnajs disks to bar con- 
figurations etc (Binney & Tremaino 1987). The most rele- 
vant analogy here is perhaps the transitions from Maclaurin 
spheroids to stationary Dedekind ellipsoids with fixed con- 
figurations in space (e.g., Chandrasekhar 1969). 

Let us first perform an analysis on aligned secular bar- 
like instability in a composite SID system. Similar to the 
single SID case of Shu et al. (2000), the criterion postu- 
lated here is expressed in terms of the ratio of the kinetic 
energy of SID rotation T to the absolute value of the gravitar 
tional potential energy W (Ostriker & Peebles 1973; Binney 
& Trcmainc 1987). In integral forms, wo derive 

T = / -T,'o{rnsf2nrdr + / -T.^oir^af'^T^rdr (65) 
Jo ^ Jo ^ 

and 

W = - r{^'o + ^'o)^2nrdr (66) 
Jo dr 

by definitions, where R is a radius that is allowed to go 

to infinity. To derive the virial theorem for the composite 

SID system from the background equilibrium conditions, we 

write equations (2) and (5) in the forms of 

-Egrn? = -^(a^Eg)-Eg^ (67) 

and 

-Sgrn^ = -A(,^Sg)-Eg^^. (68) 
Adding equations (67) and (68), we obtain 
(Sgnf + Egn^)^= |.(„2sg + a^Sg) + (Sg + Eg)^ . (69) 

Multiplying equation (69) by 2-!Tr^dr and integrating from 

to R, we arrive at 

2{T+U) + W = 27ri?'[a^ES(ii) + a^Eg(i?)] , (70) 
where 

fR 

U= (o^Eg -f OgEg)27rrdr (71) 
is the thermal energy contained in the composite SID sys- 
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tem. Equation (70) stands for the virial theorem generahzed 
for a composite SID system. 

Using cquihbrium equations (8) — (11), the two integrals 
(65) and (66) can be expressed as 

Dj + {Dl + l-l/P)6 



T = at{Di + l)- 



2G{l + 5) 



-R 



and 



atjD^s + 1)' 
G 



R 



(72) 



(73) 



For a composite SID system of infinite radial extent, both 
integrals diverge as R —^ oo but their ratio remains finite. 
So the ratio of the kinetic energy of rotation to the absolute 
value of the gravitational potential energy is 



T ^ Dj + {D'i + 1 - 1/I3)& 
|W| 2(l + Di)(l + 5) 



1 + 5/(5 



(74) 

Equations (72) — (74) can all be explicitly symmetrized with 
respect to the parameters of the two SIDs by using back- 
ground condition (11). We here use equation (74) as an ex- 
ample to illustrate this symmetry between parameters of two 
SIDs, namely, 

T ^ alT„Dl + ggSggg 

|W| [ai(l + D^) + o2(l + L>i)](S, + Eg) ■ 

Note that the value of T /\W\ falls between and 1/2 as 
usual (Binney & Tremaine 1987) and increases with the in- 
crease of DI . Therefore, for stationary aligned configurations 
in a composite system of two coupled SIDs, the two possi- 
ble values of Dl correspond to two values of T ratio; 
the larger and smaller values of Dl correspond to larger and 
smaller values of T/|VV| ratio. 

To illustrate the physical significance of the above re- 
sult, let us examine a specific case when |m,| = 2, 5 = 
0.25 and l3 = 2. Equation (51) now yields two solutions 
Dl = 1.0552 and Dj = 0.0948 (see Fig. 1). Substitution 
of these two values of Dl into equation (74) in order gives 
T/\W\ = 0.2810 and T/\W\ = 0.0890, respectively Based 
on numerical simulation experiments (involving 300— body 
particles) for stability of a rotating disk, Ostriker & Peebles 
(1973) suggested empirically that the approximate condi- 
tion T/|>V| < 0.14 ± 0.02 is necessary but not sufficient 
for stability against bar-type instabilities'. That is, when 
T/\W\ > 0.14 ± 0.02, a disk system would rapidly evolve 
into bar-type configurations (Miller et al. 1970; Hohl 1971; 
Hunter 1977). By our numerical example above and the anal- 
ogy to the single SID case (Shu et al. 2000), the transition 
criteria from an axisymmetric equilibrium to aligned secu- 
lar barlike instabilities via two different modes (i.e., upper 
yi and lower 3/2 solutions) correspond two different values of 
T/|W| ratio (larger and smaller than ~ 0.14 respectively). 
These variations of T/j W| ratio appear considerable but not 
totally surprising because for Kalnajs disks (Kalnajs 1972), 



^ Binney & Tremaine (1987) seems to suggest both necessity and 
sufficiency of this criterion. 



the relevant criterion is T/|W| < 0.1268 for the stability of 
bar modes (see also discussions of Binney & Tremaine 1987). 

There are two possible interpretations for the 1/2 solu- 
tion with a fairly low T /\W\ = 0.0890. First, if the corre- 
spondence between the ^/|W| ratio and the onset of sec- 
ular barlike instabilities holds as suggested by Shu et al. 
(2000) for a single SID case and if our extension of this cor- 
respondence to a composite SID system is valid, then our 
analysis suggests that the threshold of T/|W1 ratio can be 
lowered considerably in a composite SID system. Second, 
as emphasized by Ostriker & Peebles (1973) in their note 
added in proof, there can be other instabilities that occur 
for T /\W\ < 0.14 ± 0.02. If this is true, then the suspected 
correspondence between the T/|W| ratio and the onset of 
secular barlike instabilities might be coincidental in the sin- 
gle SID case. We are inclined towards the first interpretation 
although numerical simulations involving a composite disk 
system are deemed necessary to resolve this important issue. 
At any rate, a composite SID system tends to be stabilized 
by introducing a sufficiently massive dark-matter halo for 
both solutions (see our analysis on composite partial SID 
system later). 

For the same value of /3 = 2 but with a smaller 5 value, 
e.g., 6 = 0.1, we then have Dl = 0.8557 and Dl = 0.2125, 
corresponding to T/\W\ = 0.2428 and T/\W\ = 0.1064, 
respectively. To go further for /? = 2 and 6 = 0.05, we 
have Dl = 0.7500 and Dl = 0.2857, corresponding to 
T/\W\ = 0.2211 and T/\W\ = 0.1204, respectively By this 
sequence of three numerical examples with fixed (3 = 2 and 
m — 2, the tendency is clear that the ratio T/|VV| for yi 
decreases but remains fairly high, while the ratio T/|VV| for 
j/2 gradually approaches the usual estimate of ~ 0.14 as 5 
becomes smaller. 

Wc now complement these numerical examples by an 
analytical analysis. As defined by equation (74), the ratio of 
T/\W\ involves three parameters 3, (S and Dl explicitly, and 
increases with increasing values of Dl and /3. Once d and 
(3 are given, the marginal value of Dl can be determined 
by equation (51) for the onset of aligned barlike instabilities 
with m = 2. We investigate below variation trends of T/| W| 
versus 6 when (3 is specified. 

Let us first start with = 1.5. For m = 2, the critical 
value /3c as given by equation (56) becomes 



(3c 



1 + 



4<5-hl. 

The condition of < ft, (sec the lower branch y2 solutions 
in Fig 1) requires that for a given 0, S must fall within a 
proper range so as to render the 1/2 sohition branch physi- 
cally meaningful. It turns out that for — 1.5, S can take 
arbitrary values from to 00. We then solve equation (51) 
for aligned stationary perturbations with m = 2 and = 1.5 
to obtain the 2/2 = Dl solution branch as an explicit function 
of 5, namely 

~ 4(1 + d) 

Substitution of this y2 into equation (74) with (3 = 1.5 yields 
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Figure 2. Variation trends of the marginal ratio T/|>V| versus 
S for different values of /3 = 1.5, 2, 2.5 at a fixed value of rri = 2. 
Each curve corresponds to an allowed range of S for 2/2 = ^3 > 0. 



an expression for marginal ratio as a function of 5, 

namely 



T 



9 + - + 1(k) + Kk)-)^ 



6[7 + 85 - (1 + 165 + 165^)1/2] ' 
which has a minimum value of 5/36 « 0.1389 at 5 = 2/3 « 
0.67 (see the upper curve in Fig. 2). 

For a somewhat larger /3 such as = 2, the value of S 
is no longer arbitrary in order to obtain nonnegative 3/2; it 
turns out that < 5 < 0.5. Repeating the same procedure 
for the 13 — 1.5 case, we then derive another expression for 
critical ratio T/|W| in terms of S as 

T _ -4-115 + 3(85 + 95")^/^ 
|W| ~ 6[-4- 55 +(85 + 952)1/2] ■ 

Within the interval of < 5 < 0.5, this ratio attains a 
minimum value of 1/12 ~ 0.0833 at 5 = 0.5 (see the middle 
curve in Fig. 2). 

By similar considerations for the case oi (3 — 2.5, the 
required range of 5 becomes < 5 < 1/8, and the critical 
ratio T/|W| attains the minimum value of 1/30 at 5 = 1/8 
(see the lower curve in Fig. 2) . 

With a further increase of (3, the required range of 5 
diminishes. When /3 > 3, there is no allowed value of 5 that 
makes j/2 physically meaningful. 

These three examples are shown in Figure 2 for the vari- 
ations of the marginal ratio T/|W| versus 5 with different 
values of /? at a fixed value of m = 2. Each curve corresponds 
to an allowed range of 5 for 1/2 = > 0. Our analysis of 
aligned cases already indicate a more complicated stability 
properties of a composite SID system than those of a single 
SID system. In particular, possible modes of smallest Dg can 
have different values of m as shown in Figs. 3 and 4. 

Starting from a specific example with 5 — 0.2 and /3 = 
1.5, the j/2 solution branch of aligned equation (51) yields 
a curve of versus the aaimuthal wavenumber m > 2 as 
shown in Fig. 3. Qualitatively, this curve differs from the 
result of one-SID case, where D'^ = \rn\/{\m\ + 2) mcreases 
monotonically with increasing m > 2 (see equation (27) of 
Shu et al. 2000 and our equation (52)). Here, attains a 
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Figure 3. The solution of j/2 branch from aligned equation 
(51) versus azimuthal wavenumber m > 2 with parameters S = 
0.2 and /? = 1.5. The smallest at m = 7. 
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Figure 4. The D'^ solution of j/2 branch of from aligned equation 
(51) versus azimuthal wavenumber m > 2 with parameters S = 
0.2 and /3 = 1.2. The smallest D"^ at m = 2. 



maximum at m = 2 and 3 and a minimum at m = 7. In this 
particular case, both rn = 2 and m = 3 give the same value 
of Dg according to aligned equation (51). 

For the same 5 = 0.2 but a smaller P — 1.2, the curve 
of Dg for 2/2 solution branch of equation (51) versus m > 2 
is shown in Fig. 4 with the smallest Dg at m = 2. Note the 
variation in solution structures in Figs. 3 and 4 for a slight 
change of /3 value. 



3.2 Unaligned or spiral disturbances 

For stationary unaligned or spiral disturbances, we take the 
following set of exact density-potential relation, 

n" = cr^r"^/" exp(iQlnr) , 

9 g -3/2 /. , N C^^) 
fx =a r ' exp(ialnr) , 
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V = v^r cxp[ia\nr) + v^r exp(ia\nr) , (76) 

where cr", a®, a, u', are constant coefficients and v" and 
ct" and and cr® are related by 

= -27rGMn(a)cT'' , 

= -2nGJ\fm{a)a'' . 

Here, A/'m(a) = K{a,m) is the Kalnajs function (Kalnajs 
1971). For spiral perturbations in a composite SID system, 
this seems to be a sensible extension of earlier work rele- 
vant to the subject (e.g., Lynden-Bell & Lomos 1993; Syer 
& Tremaine 1996; Shu et al. 2000; Lou 2002). Note that the 
radial scaling parameter a (related to the radial wavenum- 
bci) is naturally taken to be the same in both stellar and 
gaseous SIDs. 

In our analysis and computations, we shall make use of 

two useful fornnila of A^,n(a). One is the recursion relation 
in m oi Nm{ci) for a fixed a (Kalnajs 1971), 



Urr.+ l{a)Urr.{a) = [(m + 1/2) V q]" 



(77) 



and the other is the asymptotic expression of Afmia) (Shu 
et al. 2000), 

J\fm{a) « (m' + a' + 1/4) "^/^ (78) 

for m? -\-o? 3> 1. When accuracy is not that crucial in some 
quantitative considerations, asymptotic expression (78) is 
also used to compute the \m\ = 1 spiral solutions later. 

For \m\ > 0, we proceed to solve equations (39) and (40) 
which, with a little algebra, can be cast into the compact 
forms of 

^2 d cP 



(79) 



d 



d' 



(^2 V + Gar-M^) , (80) 



r dr dr'^ ^ 
where coefficients Hi, H2, G\ and G-z are defined by 



Hi 

H2 = 



1 



1 



2) 



Gi - 
G2 = -- 



1 



1 - 



{Di + r).\r,n{a) 



1 + 5 

{Dl + l)Nm<,a)d 



1 + 5 



(81) 



Df(m2-2) 1 + 5 

{Dl + l) Nm(a)5 



(82) 



Dl{m?-2) 1 + 5 • 

We substitute ijf and /x® in the forms of equation (75) and 
(76) into equations (79) and (80) to obtain 

[1 - ffi(m^ +a + l/^W = Gi(m^ + a + l/4)/i« , (83) 

[1 - H2{rn +a + l/4)]/i» = G^irn +a + . (84) 

Without the gravitational coupling between the stellar and 
gaseous SIDs, the result from either disk would be the same 
as that of Shu et al. (2000) for a single SID. 

By multiplying both sides of equations (83) and (84) 



and removing /i^pi®, we obtain 

[1 - Hi{m' + a' + 1/4)] [1 - i/2(m' + a' + 1/4)] 
= GiG2(mV 0^ + 1/4)^ . 

By equations (81) and (82), equation (85) becomes 



(85) 



Df(m2 -2) 



1 + 5 



2 , 2,1 
m + a + - 
4 



D|(m2 -2) 



{Dl + l)Nm{a)5 



\ + 5 

X I + a' + \ 
4 



- ^' + 1 A/'^wy ^^ 2 1 

D2(m2-2)D|(m2-2) (1 + 5)2 V 4 

(86) 

Equation (86) can be rewritten in a physically more 

informative form of 



x\Dl{m^-2) 



_ [pl + \)M^[a) 
1 + 5 



X I + + ^ 



1 - 



(D2 + l)Mn(a)5 



1 + 5 



X I + + ^ 



(87) 



I ^^m2 , ,N AC(a)5 ( 2 , 2,1 



which is an exact relation among /3s, Dg, m, a and 5. 

By substituting relation (11), namely Z)g = fiip], + 
1) — 1, into equation (87), we obtain a quadratic equation 

oiy = Dl, 



+ Ciy + Go = , 
where three coefficients are defined by 
G2 = /?(m^ - 2)[m^ - 2 + (m^ + + 1/4)A^^(q)] , (88a) 



Gi 



f^(-^-2)(m^+a^ + i)AA„ 



x(m'+Q' + i) Um{a)-{m^ -2)(a^ + 1 



1 + 5 
/3 (886) 



2m^ + — T 
4 



^2 _ 2 + ('^^+ + V4)AA^(q) 



Go = ( + a' + J ) <! (m' - 2) 



1 + 5 
1 + 5 



x[m^ + a + \ ]Mm{a) \f3+[rn +a +^ 



1+5 
V 



x[2m^ +a - \ 



1 + 5 



(88c) 



3.2.1 The analytical case of (3 = 1 

Parallel to the analysis for aligned case, we first investigate 
the /3 = 1 case to gain useful insight. When (3 = 1 and thus 
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Dg = Dg by condition (11), equation (88) becomes 
[£),'(m^-2)-(m' + a' + l/4)] 

x{L>?[(m^ - 2) + Mri(a)(m' + + 1/4)] 
+[J\fm{a) - l](m^ + a+ 1/4)} = 
which has two branches of solution 

^ + + 1/4 



Ds=yi = — 



TV? -2 



and 



i?. = 2/2 = 



\l-Mm{a)]{m^ + 1/4) 
{m? - 2) +Mn(a)(m2 + + 1/4) 
Meanwhile from equation (83), we have 

/^s _ 1- Hi{m^ + + 1/4) 



(89) 



(90) 



(91) 



Gi(m2 + a2 + 1/4) 
[D,^(m^ - 2) - (m^ + + 1/4)] (1 + J) 



(92) 



(Df + l)Mri(a)(m2 + a2 + 1/4) 

For Dg = j/i in equation (92), one gets / = —1 and for 
Dg = j/2 in equation (92), one gets fJp / = S. These results 
parallel those of the aligned case. Note that by solution (90) , 
one has unphysical cases of < for |m| = 1. For the 
solution branch (91) with \m\ = 1, the situation is somewhat 
involved. We shall return to this case later at the end of 
subsection 5.^.5 around equations (100) - (104). 

In general, (3 > 1 and 5 varies. For different values of 
5 and /3 with \m\ > 1, equation (88) or (89) always yields 
two solutions of D^. For \m\ — 1, one solution is unphysical. 
For |m| > 2, we find that for the two solutions of Dg, the 
ratio of surface mass density perturbations jJp / jJ.^ is either 
positive or negative. This means that for stationary spiral 
disturbances in a composite SID system, fx^ and n" are either 
in phase or out of phase. We shall show details later. 

Here we have just discussed the special analytic case of 
(3=1 and point out that the \m\ = 1 case is unique. Our 
subsequent computation and analysis for general (3 values 
arc case-specific from the |?n| = case to the \m\ > 2 cases 
and then back to the Iml = 1 case. 



3.2.2 Marginal stability of axisymmetric disturbances 

While equation (88) or (89) is derived for cases of |m| > 0, 

it can also be used to describe the marginal stabihty for the 
special CEise of |m| = 0, that is, stationary cixisymmetric dis- 
turbances. Equation (89) with \m\ = yields two branches of 
solution constrained by the physical requirement" of > 0. 
For the positive portions of solution D j with not-so-extreme 
parameters, we find the basic Dg versus a profile is qualita- 
tively similar to Figure 2 of Shu et al. (2000) as expected. 
When Dl falls in the range between the maximum of the 
collapse regime and the minimum of the ring fragmentation 
curve, the composite SID system can support stable oscilla- 
tions as indicated in Figs 5 — 9. In the limit of a — > 0, stable 

II Not shown here is a second solution branch of Dg which is 
negative and thus unphysical (see Figure Dl in Appendix D for 
an example corresponding to the same parameters of Figure 2 
and also see equation (90) for /3 = 1 with m = 0). 
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Figure 5. The marginal stability curve of Dg versus a when 
m = 0, 5 = 0.2, /3 = 1.5. As S is small and /3 is slightly greater 
than 1, the curve does not differ significantly from the single SID 
case shown in Fig. 2 of Shu et al. (2000). Parameter regimes 
of collapse and ring fragmentation instabilities as well as stable 
oscillations are marked. 
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Figure 6. The marginal stability curve of Dg versus a for m = 0, 
S = 0.2, f3 = 10. While S remains small, a larger lowers the ring 

fragmentation boundary. It is easier for the system to become 
unstable in the form of ring fragmentation but with a larger a. 



"breathing modes" can be sustained (Lemos et al. 1991). 

When 5 = Sq/Sq is small, which means the gas mass is 
small compared to the stellar mass, the gas influence on the 
stellar disk is weak and the marginal stability will be slightly 
modified compared to that of a single SID situation. As 5 
increases, the role of gas disk becomes dynamically more im- 
portant. For young disk galax:ies during their early epochs 
of formation and evolution, the gas materials can be more 
abundant than or comparable to stellar mass, the axisym- 
metric marginal stability should be quite different from the 
single SID situation. 

As 5 and/or /3 vary, the marginal stability curves change 
accordingly as shown in Figures 5 — 10. For small 6 and 
f3 m 1, the shape of marginal stability curve is similar to 
that of a single stellar SID (Shu et al. 2000). The curve con- 
tains two separate branches in the Dg — a plane. The first 
collapse branch starts at the vertical Dg>0 axis with a = 
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m=0, 6 =1, e =10 




Figure 7. The marginal stability curve of versus a when 
m = 0, (5 = l,/3 = 10. The increase of S shrinks the collapse 
regime at the lower-left corner, while the increase of both S and 
/3 lowers the marginal stability curve of ring fragmentation. 
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Figure 8. The marginal stability curve of versus a when 
m = 0, 5 = 5, (8 = 3. For a moderate /3 value, the increase of S 
strongly suppresses the collapse regime while lowers the threshold 
to ring fragmentation considerably. 
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Figure 9. The marginal stability curve of versus a when 
m = 0, <5 = 10, /3 = 5. For a sufficiently large 5, the collapse regime 
disappears while the danger to ring fragmentation increases. 
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Figure 10. The marginal stability curve of versus a when 
m = 0, (5 = 1, /3 = 30. The ring fragmentation curve crosses the a 
axis (and becomes negative) such that ring fragmentation occurs 
for all when a is sufficiently large. A much reduced collapse 
regime still exists. 



and goes down to the horizontal a > axis with Dg = 0. 
The second ring fragmentation branch has a vortical asymp- 
tote, where the value of approaches infinity at a finite a. 
The curve descends with increasing a to a minimum and the 
ascends to infinity as a ^ oo. There are no other positive 
values for that can be found between the two branches 
just described. When 6 is fixed and /3 takes larger values, 
which means the stellar velocity dispersion further exceeds 
the sound speed of the gas SID, both branches descend, but 
this trend of variation is more conspicuous for the ring frag- 
mentation branch. 

For d = 0.2, the increase of /3 significantly lowers the 
ring fragmentation curve. In comparison, the collapse branch 
changes slightly with increasing /3 when S is small. This is 
shown in Figs 5 and 6. As 6 increases for fixed /3 value, 
the collapse branch shrinks with the surrounded area re- 
duced (see Figs. 6 and 7). The collapse regime continues to 
shrink with increasing 5. When 6 becomes sufficiently large, 
this branch will completely disappear from the — a di- 
agram (see Figs. 8 and 9), leaving only the ring fragmenta- 
tion branch (Fig. 9). Meanwhile for the ring fragmentation 
branch, the trend is clear that as 5 increases, the minimum 
of the marginal stability curve decreases at larger values of 
a. Eventually, the curve crosses the a axis and becomes neg- 
ative when S and /3 take proper values (Fig 10). We conclude 
that as (5 or /3 increases, both marginal stability curves move 
downward (towards smaller Dj). 

It is interesting to note that the location a of the vertical 
asymptote for the ring fragmentation branch remains fixed 
independent of 6 and /3. By equation (88), diverges when 
C2 defined by equation (88a) vanishes, or equivalently**, 
M)(a)(Q;^ + 1/4) = 2. By recursion relation (77) and the 



** This is also true when magnetic field is involved in a SID (Shu 
et al. 2000; Lou 2002; Lou & Fan 2002). 
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asymptotic expression (78), wc may choose approximately 

Mo{a) = (9/4 + q')/(1/4 + a')M2(a) 

= (9/4 + q')/[(1/4 + a")(17/4 + q')^/'] (93) 

to estimate an a value of 1.7928 for C2 = in equation (88) 
(Shu et al. 2000; Lou 2002; Lou & Fan 2002). 

Wc now provide physical interpretations for the two 
curves of marginal stability. The parameter a is a measure 
of radial wavenumber. Small a corresponds to larger radial 
perturbation scale. Therefore the composite SID system is 
vulnerable to Jeans' instability when a is sufficiently small. 
Because of the angular momentum conservation, a SID rota^ 
tion tends to work against gravitational collapse. Therefore, 
it requires a smaller a to induce Jeans' collapse in the pres- 
ence of disk rotation. When disk rotation is sufficiently fast 
{Ds larger than a critical value at a = 0), Jeans' collapse 
can be completely prevented. 

On the other hand, the composite SID system should be 
also vulnerable to axisymmetric Toomre-type instability. 
In various galactic contexts, it is of considerable interest to 
search for an effective Q parameter for a composite disk sys- 
tem (Jog & Solomon 1984a,b; Bertin & Romeo 1988; Kenni- 
cutt 1989; Romeo 1992; Wang & Silk 1994; Elmegreen 1995; 
Jog 1996; Lou & Fan 1998b). Shu et al. (2000) noticed that 
the minimum of the ring fragmentation branch is effectively 
related to the Toomre Q parameter (Toomre 1964) in a very 
close manner. In a magnetized SID with a coplanar mag- 
netic field, Lou (2002) and Lou & Fan (2002) noticed that 
that the minimum of the MHD ring fragmentation branch 
is effectively related to the generalized MHD Qm parameter 
(Lou & Fan 1998a) in a very close manner. Physically, we 
therefore expect that the minimum of the ring fragmenta- 
tion branch should be intimately related to an effective Q 
parameter in a composite disk system (Elmegreen 1995; Jog 
1996; Lou & Fan 1998b). When exceeds the minimum of 
the ring fragmentation branch, the composite system of two 
coupled SIDs becomes unstable to ringlike fragmentation for 
the range of radial wavenumbers a between the two points 
of intersection of the horizontal line and the the marginal 
stability curve. 

3.2.3 Stationary logarithmic spiral configurations 

Following the same procedure of obtaining solutions for ax- 
isymmetric disturbances, we derive solution curves for cases 
of |m| = 1,2,3.... As expected, for each set of {m, 5, /?}, 
there are two solution curves for Dg versus a in a compos- 
ite SID system according to equation (88) or (89). In some 
cases, part of the lower solution may become negative which 
is unphysical. Figures 11 — 13 show a set of numerical ex- 
amples with |m| = 2. 

In the case of |m| = 1, only one branch of is physi- 

tt The original criterion was derived for local axisymmetric per- 
turbations in a rotating disk (Safronov 1960; Toomre 1964). 



in=2, 6 =0. 2, 3 =1. 5 




Figure 11. Two solutions of stationary spiral configurations for 
versus radial wavenumber a when \m\ = 2, S = 0.2, /3 = 1.5 
as derived from equation (88). For small 5 and /3 slightly greater 
than 1, the lower branch is almost the same as Fig. 3 of Shu et 
al. (2000). The upper branch is novel due to a composite disk 
system. Both curves increase with increasing a. 
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a 

Figure 12. Two solutions of stationary spiral configurations for 
versus radial wavenumber a when \m\ = 2, 5 = 0.2, /3 = 3. 
As /3 increases, both branches move downward. It is possible for 
the lower branch go across the horizontal a axis, while the upper 
branch will always remain positive. 

cally meaningful. We will discuss more on this problem later. 
Let us first examine cases of |m| > 2. 

For given values of S and /3, there are two solution curves 
in general for versus a when \m\ > 2. The lower branch 
may cross the horizontal a axis when S and/or /3 are suffi- 
ciently large. In contrast, the upper branch will never be- 
come negative. We note that variations of /3 and S play 
distinctly different roles. The increase of S mainly moves 
the upper branch upward and the lower branch downward, 
whereas the increase of p tends to move both branches down- 
ward. For very small S values, different P values may lead to 
quite different solution configurations (see Figs. 11 — 13). 

In order to fully understand solution behaviors (Figs. 
11 — 13), we perform analytic analysis on properties of 
stationary spiral configurations as in the aligned cases of 
\m\ > 2. Firstly, through extensive numerical experiments. 
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Figure 13. Two solutions of stationary spiral configurations for 
versus radial wavenumber a when |m| = 2, (5 = 1, /3 = 3. 
By increasing 5 at fixed p (see Fig. 12), the lower branch moves 
downward while the upper branch moves upward. 



one can show that y = increases with wavenumber a for 
both upper and lower branches. Secondly, once |m| and 5 
are known with fixed a, we would like to know how y = 
and / vary with /?. It turns out that for both solution 
branches, y = monotonically decrease with increasing [3 
and decreases with increasing D^. It is then clear that 

l-i^ /n" goes with f3 in the same trend. In more informative 
forms, we derive two useful inequalities for yi and j/2 below. 
The first one for yi is thus 

(m^+a^ + l/4)/(m^-2) 

[m^ - 2 + (m^ + + i/4)]A/-„(q) 



K - 2 + (m2 + a2 + l/4)J\fm (a)] ( 1 + 5) .' ^^^^ 

< yi < (m'' + l/4)/(m^ - 2) 

with the lower bound on the left-hand side determined by 
taking f3 —>■ oo and the upper bound on the right-hand side 
obtained by taking /? = 1; and the second one for j/2 is 

[l-A/'„(Q)](m^-f Q"-hl/4) 



-1< 1/2 < 



(95) 



- 2 + J\fm{a){m'^ + a"^ + 1/4) 
with the lower bound on the left-hand side determined by 
taking oo and the upper bound on the right-hand side 
obtained by taking (3 = 1. One can readily show that the 
lower bound of j/i is positive definite as Nm{a) < 1 for 
\m\ > 2. Therefore, the = j/i branch is always physically 
valid given values of S and (3 with \m\ > 2. Moreover, as in 
the aligned cases, the left-hand side of inequality (94) for j/i 
remains always larger than the right-hand side of inequality 
(95) for j/2 such that the two solution branches yi and j/2 
will never intersect with each other because j/i > 3/2 always. 
Meanwhile, there exists also a critical value (3c for y2 = 
branch at a given wavenumber a such that 2/2 = becomes 
zero. In analytic form, this critical value (3c is given by 
™2 2 + (m^ +a^ + 1/4) 



2-FAA„(a)(m2 + Q2 + l/4) 



H- 



AC(a)[l -A/"m(a)](m' -^g^ + 1/4) 
Mn(a)(m2 + a-' + 1/4)6 + {'m? - 2)[1 + 5- Mm{a)] 



Since Mm,(o.) < 1 for |m| > 2, the critical value (3c remains 
always greater than 1 and decreases with increasing 5. When 
5 — > 00, the critical value /3c approaches a limit of 

^ - 2 + (m" + + 1/4) 

PoLim - ^2 _ 2+Nm{a){m? + a^ + 1/4) " ^^'^ 

In parallel to the study of aligned cases, we further derive 

the range of /i^/zi" for both solution branches when \m\ and 
5 are known at a given wavenumber a. For the t/i = Dl 
branch, we thus have 

-{rn' + a' + l/4)N^{a)S 



-1 < 



/x" (m2-ha2 + l/4)AC(a)5-|-(m2-2)(l-F5) 



(98) 

with the lower bound on the left-hand side determined by 
(3=1 and the upper bound on the right-hand side deter- 
mined by /3 — » 00; for the j/2 = solution branch, we have 



5 < — < 



1 + 5 



(99) 



Nm{a) 

with the lower bound on the left-hand side determined by 
(3=1 and the upper bound on the right-hand side deter- 
mined by the critical value (3c such that y2 = = 0. Since 
/n" increases with increasing (3, we always have negative 
l-jp / if for the 2/1 = DI branch and positive jJp / fjC for the 
2/2 = Dt branch. 

It is remarkable that the spiral cases are strikingly 
similar to the aligned cases. If one replaces Mrn(a) by 
(m^ + + l/A)-^''^ and (m^ + + 1/4)"^/^ by all 
the above results for the spiral cases will degenerate to the 
aligned cases. By these considerations of aligned and spiral 
cases with \rn\ > 2, this apparent similarity should not be 
surprising because the common physical nature of station- 
ary aligned and spiral disturbances, that is, purely azimuthal 
and general density waves in a composite SID system bal- 
anced by disk rotation. 

As promised earlier, we now get back to the \m\ = 1 
case which is special for a stationary spiral configuration as 
well as for a stationary aligned configuration. For aligned 
\m\ = 1 case, equation (51) is satisfied for arbitrary > 
(see Shu et al. 2000 for a single SID). For spiral \rn\ = 1 case, 
when 6 and a are given, the j/i and 1/2 solution branches will 
increase and decrease with increasing (3, respectively. When 
(3 —>■ 00, the two solutions are 



y = -l 



and 



2 , 5\ A/-i(a)(a^ + 5/4)(g^ + l/4) 
" ^ 4) ^ [J^i{a){a^ + 5/4)-l]{l + 5) ■ 

(100) 

Note in the above equation, we have J\fi{a){a'^ + 5/4) — 1 
always positive according to equations (77) and (78). 

The determination of the specific bounds of these two 
solutions will depend on a relation between 5 and a. When 

2 5\ Afi{a)ia^ + 5/4){a' + 1/4) 
^ 4) ^ [Afi{a){a^ + 5/4)- 1]{1+ 5) 



(96) 



/ 2 , n <?^(<.i)(<.r+5/i)-i]-l 
\ ^ 4j lJ^i{a){a^ + 5/4)-l]{l+S) 



> -1 
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which means d < [J\fi{a){a^ + 5/4) — 1] ^ , we have 

-(a' + 5/4) <yi < -1 (101) 

and 



2 , 5\ A/-i(q)(q" + 5/4)(q^ + 1/4) 



<y2< 



[l-M(Q)](a' + 5/4) 



(102) 



M(a)(a2 + 5/4)-l ' 
On the other hand, when 

/ , 5\ M(a)(a^ + 5/4)(a2 + l/4) 
y V 4y'^[M(a)(a2+5/4)-l](l + 5) 

1 



= -1 



2 1\ 3[M(a)(a^ + 5/4)-l] , 
4; [M(Q)(a2 + 5/4) - 1](1 + 5) ^ 



which means S > [A/i(a)(a'^ + 5/4) — 1] ^ , we have 



2 5 2 
-a - - < yi < -a 



5 M(Q)(a^ + 5/4)(a^ + 1/4) 
4"^ [M(a)(a2 + 5/4) + 5) 

(103) 



and 



-1< t/2 < 



[l-A/-i(a)](a2 + 5/4) 



(104) 



M (Q)(a2 + 5/4) - 1 ■ 
Note that for bounds of these inequahties, one side is de- 
termined by /3 = 1 and the other side by /? — > 00 for both 
solution branches. It becomes clear that spiral \m\ — 1 case 
differs from spiral |m| > 2 cases. In particular, the ranges of 
yi and j/2 solutions are controlled by a relation between 5 
and Q, that is, whether 5 > [N'i{a)(a'^ + 5/4) — 1]~^ or not, 
such that the two solution branches do not intersect with 
each other. By inequalities (101) and (103), the yi solution 
branch is always negative and thus unphysical. In contrast, 
by inequalities (102) and (104), it is still possible for por- 
tions of the j/2 solution branch to be positive when 5, a and 
f3 take on proper values. Note that this j/2 solution branch 
is the counterpart of a single SID case. 

As indicated in Figs. 11 — 13, for given values of |m|, 
P and S, it possible for a composite SID system to sustain 
upper (smaller a) and lower (larger a) solution branches 
with different values of a simultaneously. Likewise, one may 
switch the role of |m| and a, that is, given a, a composite 
SID system can support various spiral solutions with differ- 
ent values of \m\ (see Fig. 3 of Shu et al. 2000). This is to be 
understood that both \m\ and a are characteristic proper- 
ties of density waves rather than background parameters of 
a composite SID system. In the linear regime, this allowed 
superposition of different solutions can give rise to a variety 
of galactic pattern structures that cannot fit by one simple 
spiral pattern. 

The interpretation for the stability properties of 
these stationary logarithmic spiral configurations is not so 
straightforward even for a single SID case because of the 
power-law background divergence at r ^ and thus the dif- 
ficulty of imposing proper boundary conditions there (Zang 
1976; Toomre 1977; Lynden-Bell & Lemos 1999; Goodman 
& Evans 1999; Shu ct al. 2000; GaUi et al. 2001). Shu et al. 
(2000) speculated that the stationary spiral solution condi- 
tion represents onset of spiral instabilities on the ground that 
transmission and overreflection of leading/trailing spiral 



density waves across the corotation can be swing-amplified. 
This was done mainly through numerical experiments and 
appeared to be consistent with an earlier proposed criterion 
(Goodman & Evans 1999). Furthermore, the stability of sta- 
tionary spiral configurations may be examined in light of the 
Ostrikcr- Peebles criterion (Shu ot al. 2000). It is then plausi- 
ble to extend this interpretation to a composite SID system 
where apparently more leading/trailing spiral density waves 
across the corotation would be involved in the processes of 
transmission and overreflection. 

Therefore, in Figs. 11 — 13, as increases from and 
first intersects with the lower solution branch, the lower 
branch spiral solution can be sustained and be vulnerable 
to dynamical barred-spiral instabilities. As increases fur- 
ther to intersect the upper solution branch, the upper branch 
spiral solution (together with the lower branch solution at 
a larger a) can be sustained and be vulnerable to spiral in- 
stabilities. In the case of Fig. 11, may be small enough 
to avoid dynamical barred-spiral instabilities. But in cases 
of Figs. 12 and 13, dynamical baxred-spiral instabilities just 
cannot be avoided for > when a is sufficiently large. 

For the specific case of Fig. 11 with m = 2, 5 = 0.2 and 
(3 = 1.5, the two minima of the yi and 3/2 solution branches, 
i.e. the values of D'^ when a ^ 0, arc 1.3406 and 0.3697, 
respectively. The corresponding ratio of T/|>V| defined by 
equation (74) are 0.2982 and 0.1552 in order. In analogy to 
the single SID case of Shu et al. (2000), the two minima 
of signal onsets of dynamical barred-spiral instability 
in the composite SID system. By this interpretation, when 
the rotation of the stellar SID _Df is smaller than 0.3697, 
the composite SID system can be stable against m = 2 spi- 
ral instabilities. Otherwise, dynamical barred-spiral insta- 
bilities will develop in the composite SID system, similar to 
processes shown in sections 5.7 and 7 of Shu et al. (2000). 
Notice the introduction of the gaseous disk has decreased 
somewhat the marginal value of Z)J and thus the threshold 
ratio of T/|Wj in comparison to a single SID case. 

Similar to the aligned cases, spiral stability properties 
of a composite SID system are more complicated than those 
of a single SID. For stationary logarithmic spiral configura- 
tions, Dg usually attains a minimum value at a = for both 
upper and lower solution branches of spiral equation (88), 
except for cases when 2/2 = Dg solutions become negative 
for small a values (see Figs. 5-9). 

As in the aligned case, wo start from the same parame- 
ter set of (5 = 0.2 and P = 1.5 yet with a fixed a = 0. The j/2 
solution branch of spiral equation (88) then yields a curve 
of Dg versus m > 2 as shown in Fig. 14 where the smallest 
value of = y2 occurs at azimuthal periodicity m = 9. 

For the same 5 = 0.2 and a = yet a smaller (3 — 1.2, 
the curve of Di for 2/2 solution branch of spiral equation (88) 
versus m > 2 is shown in Fig. 15, where m = 2 gives the 
smallest D^. 
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Figure 14. The solution of 2/2 brancli from spiral equation 
(88) versus azimuthal wavenumber m > 2 with parameters a = 0, 
S = 0.2 and P = 1.5. The smallest D'^ at m = 9. 
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Figure 15. The solution of j/2 branch from spiral equation 
(88) versus m > 2 with a = 0, S = 0.2 and /3 = 1.2. The smallest 
D? at m = 2. 



4 COMPOSITE PARTIAL SID SYSTEM 

In this section, we consider stationary perturbations in a 
composite system of two gravitationally coupled partial SIDs 
(Syer & Tremaine 1996; Shu et al. 2000; Lou 2002; Lou & 
Fan 2002) to model the gravitational influence of a back- 
ground axisymmetric dark-matter halo on the composite 
SID system. The response of this massive dark-matter halo 
to perturbations in the composite SID system is ignored. 
Formally, we introduce am additive gravity term d^/dr or 
d^/dip to the basic equations (2), (3) and (5), (6), where $ 
represents the gravitational potential contribution from the 
background dark-matter halo. 

For momentum equations in the stellar disk, we have 



at 



du" f dw 



dr 



dip 



1 9 2 



0^ _ a£ 

dr dr 
(105) 



and 

Likewise, for momentum equations in the gas disk, we have 



dt 

and 

dt 



dr 



3 



92 



dji 

dr 



dip 



dip 



'S9 dr 



(a^E«). 



dr 



_ 9£ 
dr 
(107) 



1 9 I 2^9^ 



For a composite partial SID system, we now introduce a 
dimensionless T parameter defined hy = 0/(0 -|- $). Now 
the background equilibrium should bo modified accordingly. 



As before, 
sities now become 



we write = agDa/r, f2g = agDg/r, and thus 



The equilibrium surface mass den- 



•y9 _ 

^0 — 



a.. 



2-KGr 

' (1 



a„ 



1 + 5 

- Dl)T5 



(109) 



2TiGr 1 + 5 

Note that condition (11), namely, al{Dl -|- 1) = dl{Dl + 1) 
still holds. We refer to partial SIDs for < < 1, in contrast 
to fufl SIDs with T ^1. 

The linearized coplanar perturbation equations should 
take the same forms of those in full SIDs written down in 
section 2.3. We thus follow the procedure of the full SIDs 
case up to equations (37) and (38). For the modified equi- 
librium with w = 0, we have for the stellar partial SID 



M + 



Dl{m^ - 2) 



^[rii + 



dr 



dr' 



TV 



27rG \ + 5 



(110) 



and for the gaseous SID 
1 



D2(m2 



r/i" -I- 



dr 



dr^ 



1 + Dg 5rv \ 



(111) 







27rG 1 + 5) _ 

(see equations (39) and (40)). Equations (110) and (111) are 
to be solved simultaneously with Poisson's integral (32). 

We now consider consequences of a composite partial 
SID system in comparison with a composite system of two 
full SIDs. 

For aligned cases: the \m\ = disturbance is still a 

rescaling of one axisymmetric equilibrium to a neighbor- 
ing axisymmetric equilibrium. For aligned cases of \m\ = 1, 
the generalized counterpart of equation (51) yields only one 
degenerate solution y = —1 which is unphysical. In other 
words, eccentric displacements of |m| = 1 are not allowed 
in a composite partial SID system. With 5 and |m| given, 
\m\ > 2 cases yield two solutions satisfying the following 
inequalities 

2\m\{r,i^ - 1)T 



2 (m2-2)(m2-h|m|J^-2)(l + <5) 
< yi < nx' l{ra' — 2) 



(112) 
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and 

-1< j/2 < |m|(|m| - J^)/(m" + \m\T - 2) , (113) 

respectively (see inequalities (54) and (55) for the cases of 
full SIDs) . For a composite system of two coupled full SIDs 
with = 1, these results naturally reduce to the familiar 
ones, namely, inequalities (54) and (55), derived earlier. 

Given parameters 5 and /3 with \rn\ > 2, partial SIDs 
lead to two solution branches of that are larger than 
those of full SIDs. Therefore, the critical value /3c for j/2 = 
= is larger in partial SIDs, namely 

l)/(m^ + |m|JP"- 2) 

|m|(|m|-l)^ 



/3c = 2(m^ 
1 + - 



|m|(m2 + \m\J^ - 2)6 + (m^ - 2)(|m| - JF) 

(114) 

For full SIDs with = 1, expression (114) reduces to ex- 
pression (56) as expected. Meanwhile, the ratio of surface 
mass density perturbations /j,^ /fi" becomes 

!^^_-^_ [DUm^-2)-m^]{l + S) 



|m|(/)| + l)^ 

Specifically for the two solution branches of Dg. 
two inequalities, namely 



(115) 
we derive 



-1 < 



< 



(116) 



(117) 



II' (m2 + |m|JF-2)(5 + m2-2 
for the = j/i solution branch, and 

5< ^ < \m\{l + 6)/T-l 

for the DI = branch. Notice that expressions (115), (116) 
and (117) have their obvious counterparts of the full SIDs 
cases of expressions (53), (58) and (59) derived earlier. 

The same procedure can be applied to the unaligned 
situation, which is expected to yield the same trends as the 
aligned situation. 

We now derive the virial theorem in a composite partial 
SID system. By equations (105) and (107) for the equilib- 
rium state, we obtain 

d{<j>o + $) 



dr 



(isSo) 



and 



y9 
^0 



dr 



d{4>o + $) 



(118) 



(119) 



dr'^ " dr 

Adding equations (118) and (119), multiplying the resulting 
equation by 2-Kr^dr and integrating from to R, we obtain 

2{T + U)+W = 2'KE?[alY.t,{R) + alT.l{R)] , (120) 

where 



(rQs) 2'Krdr + 



isg(rng)^27rrdr , (121) 



rR 

U= (o^Eg -|-OgEg)27rrdr , (122) 
Jo 



■r.(Eg + Eg)^(^2..d. 



dr 



(123) 



are the kinetic energy of rotation, the thermal energy and 
the gravitational potential energy, respectively. Note here 



the dark-matter halo contributes to the gravitational poten- 
tial energy and = 0o/(0o + *)• 



5 SUMMARY AND DISCUSSIONS 

There are four main tracks of research that lead to our inves- 
tigation here on stationary perturbation configurations in a 
composite system of two gravitationally coupled full or par- 
tial SIDs. The first track is the classical study on spheroidal 
and ellipsoidal figures of equilibrium (Charidrasekliar 1969; 
Tassoul 1978) named after famous mathematicians includ- 
ing Maclaurin, Jacobi, Dedekind, Riemann and Poincare. 
These spheroids and ellipsoids have their analogous com- 
pressible counterparts (Ostriker 1978), disk counterparts by 
collapsing along the rotation axis (Weinberg & Tremaine 
1983; Weinberg 1983), as well as counterparts of composite 
systems (Vandervoort 1991a,b). The second track is along 
the past investigation of possible SID configurations and 
their stability properties (Zang 1976; Toomre 1977; Lemos, 
Kalnajs & Lynden-Bell 1991; Lynden-Bell & Lemos 1993; 
Lee & Goodman 1999; Evans & Read 1999; Goodman & 
Evans 1999; Syer & Tremaine 1996; Shu et al. 2000; GaUi et 
al. 2001). The third track follows the studies of density waves 
in a composite disk system (Lin & Shu 1966; Jog & Solomon 
1984a,b; Bertin & Romeo 1988; Elmegreen 1995; Jog 1996; 
Lou & Fan 1998b, 2000a,b). The fourth track pursues a ba- 
sic understanding of MHD density waves in magnetized disk 
systems (Fan & Lou 1996; Lou & Fan 1998a, 2002, 2003; Shu 
et al. 2000; Lou, Yuan & Fan 2001; Lou 2002). Naturally, 
we shall further study possible stationary configurations in 
a composite system of gravitationally coupled stellar SID 
with gaseous MSID. These studies are of interests in galac- 
tic dynamics (Bertin & Lin 1996), star formation (Shu et al. 
1999), the light cusps seen in the nuclei of galax;ies (Crane et 
al. 1993), circumnuclear starburst "rings" and disk accretion 
processes axound active galactic nuclei (Lou et al. 2001). 

In reference to the work of Lou & Fan (1998b) on den- 
sity waves in a composite disk system, of Shu et al. (2000) 
on stationary perturbation configurations of a single isope- 
dically magnetized SID, and of Lou (2002) on stationary log- 
arithmic fast and slow configurations in a single MSID, we 
studied here possible stationary configurations in a compos- 
ite system of two gravitationally coupled full or partial SIDs 
using the two-fluid formalism, solved the Poisson equation 
exactly, and derived both aligned and unaligned or spiral 
stationary perturbation solutions. We have reached several 
conclusions and suggestions summarized below. 

(i) Aligned Stationary Configurations 

For perturbation configurations in a composite SID sys- 
tem, there exist in general two branches of solutions as ex- 
pected (Lou & Fan 1998b). This is one major difference be- 
tween a composite SID system and a single SID case. We 
shall start with the aligned cases first. 

The aligned axisymmetric |m| = case is again a rescaling 
from one axisymmetric equilibrium to a neighboring axisym- 
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metric equilibrium (Shu et al. 2000; Lou 2002), except that 
now the rescaUng occurs simultaneously in both SIDs. 

Wc find as in the single full SID case (Shu ct al. 2000) 
that, for eccentric |m| = 1 displacements, equation (51) can 
be satisfied for unconstrained positive values of D^. Given 
known idealizations of the full SID model, wo consider that 
such eccentric displacements are possible. However, the situ- 
ation is quite different for a composite system of two partial 
SIDs for which eccentric displacements are not allowed con- 
sistent with the similar result of a single partial SID (Lou 
2002). 

For |m| > 2, we derive two branches of solution for 
possible values of rotation speed parameter such that 
aligned stationary configurations are sustained in a compos- 
ite SID system. Of the two branches, one is always larger 
than the other. For the larger branch denoted by yi , pertur- 
bation enhancements of surface mass densities in stellar and 
gaseous disks are out of phase, while for the smaller branch 
denoted by j/2, the two enhancements are in phase. By the 
physical requirement of > 0, the larger j/i branch is al- 
ways valid (see inequality 54), while the smaller y2 branch 
may become negative and thus unphysical (see inequality 
55). Given \m\ > 2 and 5, when /3 exceeds a critical value 
/3c defined by equation (56), the smaller j/2 branch turns 
negative (see Fig. 1). In other words, only for disk galaxies 
with /? < /3c, the smaller j/2 branch can be physically mean- 
ingful. Moreover, the smaller y2 branch corresponds to the 
solution of single SID case as it would approach the limit of 
the single SID result when /3 — » 1, that is, when the system 
can be effectively viewed as a single SID. In contrast, the 
larger j/i branch obtained in a composite SID system has 
no counterpart in the single SID case. Physically, stationary 
aligned disturbances actually involve purely azimutlial prop- 
agation of density waves counter balanced by SID rotation 
(Lou 2002; Lou & Fan 2002). The same density- wave inter- 
pretation holds for stationary unaligned or spiral perturba- 
tions. As in the single SID case, stationary aligned configu- 
rations here mark onsets of secular instabilities from an ax- 
isymmetric equilibrium to nonaxisymmetric configurations. 
For different values of |m|, distinct aligned configurations 
may coexist. 

(ii) Unaligned or Spiral Stationary Configurations 
Solution properties of the unaligned or spiral logarith- 
mic configurations bear strong resemblance to those of the 
aligned cases. The spiral case involves both radial and az- 
imuthal density wave propagations, whereas the aligned case 
involves only an azimuthal propagation of density waves. 
To a certain extent, wc may regard aligned perturbation 
cases as special examples of spiral perturbation cases with 
\m\ > 2. 

For the marginal case of \m\ = 0, we have computed the 
marginal stability curves of versus radial wavenumber 
a, invoking the same physical interpretations of Shu et al. 
(2000) for regimes of collapse, stable oscillations, and ring 
fragmentation. We shall further demonstrate the correctness 
of this interpretation by performing a WKBJ analysis in a 
separate paper. For different parameters 5 and /3, we find 



these marginal stability curves may vary significantly in ref- 
erence to those obtained by Shu et al. (2000). Numerical 
examples arc illustrated in Figs. 5—10 (sec also Fig Dl 
in Appendix D for a typical solution structure). Generally 
speaking, the introduction of a gaseous SID tends to de- 
crease the stability of the composite SID system. Meanwhile, 
the possibility for the system to collapse may be suppressed 
(Fig. 9) . In primordial disk galaxies where gas disk mass ex- 
ceeds stellar disk mass, the possibility of large-scale collapse 
may be suppressed while the chance of ring fragmentation 
increases. 

The \m\ — 1 case is again similar to the single SID prob- 
lem. Specifically, the j/i solution branch, which is the novel 
solution of a two-SID system, remains always negative and 
is thus unphysical. Meanwhile the 1/2 solution branch, which 
is the counterpart of one-SID problem, can still have phys- 
ically valid portions constrained by inequalities (102) and 
(104). 

For \m\ > 2 cases, we have obtained analytical results 
almost in the same forms of aligned cases. This can be 
clearly seen by comparing equations (94) — (99) with equa- 
tions (54) — (59). By replacing the Kalnajs function A/'m(a) 
with (m^ +a^ + l/4)-i/2 and (m^ + a'' + l/4)i/2 with \m\ 
in results of spiral cases, the resulting equations become the 
same as the aligned cases. To reiterate, the aligned case is 
just a special one complementary to spiral cases, all of which 
involve hydrodynamic density waves in a composite system 
of two coupled SIDs. 

It should be mentioned that in the WKBJ regime, such 
stationary logarithmic spiral density wave perturbations ac- 
tually carry an outgoing angular momentum flux which is 
constant in r. With our composite SID model, a source of an- 
gular momentum is thus required at the origin r — where 
the surface mass density diverges as ~ r^^. More gener- 
ally, the outward advective angular momentum flux (stellar 
and gas SIDs together) is constant in r, while the gravity 
torque term involves a nontrivial integral in z (see Section 
4 of Fan & Lou 1999). In various contexts, several authors, 
e.g. Lynden-Bell & Kalnajs (1972), Goldreich & Tremaine 
(1978), and Fan & Lou (1999), have previously investigated 
angular momentum fluxes carried by density waves or MHD 
density waves in the tight-winding or WKB.1 regime. 

(iii) A Composite System of Two Partial SIDs 

In contexts of disk galaxies, the gravitational effect of the 
more massive dark-matter halo may be incorporated using 
a model of so-called partial SID by introducing a parame- 
ter which is the ratio of disk potential to the total potential 
(0 < < 1). The background surface mass densities of 
the two coupled SIDs are modified. For stationary pertur- 
bations in a composite partial SID system (0 < ^ < 1), 
it is straightforward to make use of the full SID formalism 
(.F= 1) by replacing 1/(1 + 5) with J^/{l + 5) and 5/(1 + 6) 
with J^5 /{I + 6) in equations (39) and (40). 

Naturally, all the previous full SID results are modified by 
the parameter quantitatively. In particular, the aligned 
|m| = 1 case is no longer a solution for arbitrary positive 
D^, that is, such stationary eccentric displacements are for- 
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bidden in a composite system of two partial SIDs. Also, the 
critical value (3c for physically valid j/2 solution branch (both 
aligned or unaligned) with |m| > 2 is increased somewhat. 

While our model formulation is highly idealized in sev- 
eral aspects, it is one step closer to a realistic disk galaxy. 
The theoretical results here provide a few useful concepts 
and possible clues for diagnostics of disk galax;ies composed 
of stars and gas coupled by mutual gravity on large scales. 
Analytical solutions derived here offer valuable benchmarks 
for numerical codes designed for N-body computations with 
fluid gas dynamics. The analysis here also paves the way to 
study theoretical model problems of stationary configura- 
tions in a composite system of a stellar SID gravitationally 
coupled with a gaseous MSID (Lou 2002). It is of consider- 
able interest to extend the current analysis into the nonlinear 
realm (Galli et al. 2001) to explore the large-scale dynamics 
of spiral shocks. 

(iv) Stability of a Composite SID System 

For both aligned and spiral instabilities, a composite SID 
system is more complicated than a single SID system in sev- 
eral aspects. Firstly, roughly speaking, the number of pos- 
sible stationary solutions is doubled except for unphysical 
j/2 solutions that must be excluded. Secondly, there are two 
more dimensionless parameters /? and S that may cause vari- 
ations in solution structures. Thirdly, because a composite 
SID system tends to be less stable, solutions of lower 
y2 branches may correspond to T/|W| ratios considerably 
less than ~ 0.14 (Ostriker & Peebles 1973). Fourthly, in cer- 
tain parameter regimes of 5 and (3 with m = 0, the collapse 
regime of spiral cases (Shu et al. 2000) may disappear. We 
shall further examine in detail axisymmetric stability prop- 
erties of a composite SID system in a separate paper. Fifthly, 
because of variations in solution structures, one needs to be 
very careful to determine the most vulnerable configurations 
with smallest values. We propose to test the Ostriker- 
Peebles criterion for a composite disk system by numeri- 
cal simulation experiments and expect considerably lower 
T/|VV| ratio than ~ 0.14 in certain parameter regimes stud- 
ied here. 

(v) Dynamics of Magnetized SIDs 

For the evolution of an individual spiral galaxy, star for- 
mation is one major process (e.g., Kennicutt 1989) and the 
global star formation rate is a key parameter. In the density 
wave scenario (Bertin & Lin 1996), spiral arms of higher gas 
concentration are locations of active star forming regions on 
smaller scales. Therefore, a protodisk galaxy begins with a 
gas disk will eventually evolve into a stellar disk. In between 
the two extremes, the disk galaxy involves a gas disk and a 
stellar disk with a variable mass ratio. As time goes on, their 
mass ratio as well as the ratio of velocity dispersions change. 
At various stages, large-scale instabilities may occur in a 
composite disk system to alter the evolution track. Further- 
more, the presence of magnetic field in the gas disk may give 
rise to magnetohydrodynamic (MHD) density waves (Fan & 
Lou 1996; Lou & Fan 1998) and instabilities that can signif- 
icantly affect the star formation rate. 

Magnetic field effects are not included at present merely 



for simplicity. In disk galaxies, the gaseous disk is magne- 
tized with energy densities of magnetic field, thermal gas and 
cosmic rays being comparable (Lou & Fan 2003). We derive 
such information based on synchrotron radio emissions from 
spiral galaxies. Large-scale magnetic fields typically lie in the 
disk plane, but there arc cases with out-of-planc magnetic 
fields. On various scales, MHD will play an important role in 
the gaseous disk and provide valuable dynamical as well as 
diagnostic information (Fan & Lou 1996; Lou & Fan 1998a; 
Lou et al. 2001, 2002; Lou 2002). In the current contexts, 
we simply point out two possible extensions of this theo- 
retical analysis. One is to consider stationary perturbation 
configurations in a composite partial SID system in which 
the gas disk is isopedically magnetized (Shu et al. 2000). It 
is then possible to study the influence of magnetic field on 
the stability of the system. Such a problem is potentially 
important to the development of galactic wind models. The 
other is to consider stationary perturbation configurations 
in a composite partial SID system in which the gas disk is 
magnetized with coplanar magnetic fields (Lou 2002; Lou & 
Fan 2002). The presence of magnetic fields allows for more 
dynamic freedom in the gas disk and gives rise to fast and 
slow MHD density waves (Fan & Lou 1996; Lou et al. 2001). 
We shall pursue these more realistic yet challenging MHD 
model problems in separate papers. 
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APPENDIX A: 

For the convenience of reference, we here outHne a general 
proof that equation (51) for stationary aligned perturbations 
always have two branches of real solutions of y = when 
|m| > 2. Note that m = is a rescaling of the background 
and m = 1 leads to unconstrained D^. By multiplying a 
factor (1 + 5), the coefficients of quadratic equation (51) can 
be rewritten in the following forms: 

Ca =/?(m2-2)(m-l)(m + 2)(l + 5) , (Al) 

Ci^-2{m^l)[{m + 2)3-{m^ -2)]I3 

- 2(m^ - l)[(m^ - 2)(5 + + m - 2] , 

Co = -m(m - l)[m(m + 2)S + - 2]f3 

+ 2m(m — l){mS -f m — 1) , 

where the common factor (m — 1) is made explicit in all 
three coefficients. The determinant A of equation (51) is 
then given by 



A 



Ci — 4C2C0 



= 4(m'--l)'(^/3^ + S/3 + C) , (A4) 
where three coefficients A, B, and C are defined by 



^ = [m - 2 + (m^ + m - 2)5f 



B = - 2{rn - 2){rn +m- 2){5^ + 1) 
- 2(2m'' + 2m^ - Qrn - 4m + 8)5 



C = [m^ + m-2 + {m^ - 2)5f . 
We can now rewrite equation (A4) in the form of 



A = 4(m^ - 1)^ 



4:A 



(A5) 



(A6) 



(A7) 



(A8) 



Since both ^ > and 



Am^{m^ - 2){ra^ -f m - 2)5(1 + Sf 



> (A9) 



[m? -2 + {m? +m- 2)Sf 
for m > 2, the determinant of equation (51) is clearly always 
positive. Therefore, equation (51) will give two branches of 
real solutions for _DJ = y for m > 2. 

In the same spirit, we can show numerically that for 
unaligned or spiral cases of \m\ > 1, equation (88) always 
has two branches of real solutions for = y. To avoid rep- 
etition and redundancy, we omit that demonstration here. 



APPENDIX B: 

We demonstrate here that when \m\ > 2, the two branches of 
solution of equations (51) and (88) monotonically decrease 
with increasing P as shown in Fig. 1. 
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Let us consider equation (51) for aligned cases with 
fixed |m| to simplify the computation. For example, with 
\m\ = 2 in equation (51), the three coefficients of equation 
(51) as modified in Appendix A become 

C'^ = 2(1 + 5)13 , (Bl) 



(B2) 



C'^ = -{4S + l)f3 + 3{2S + l) (B3) 

by further removing a common factor 4. Explicitly, the two 
solutions of equation (51) can then be written as 



2/1,2 = 



2C- 



(B4) 



Substituting expressions (Bl) - (53) for C2 , C" and d',' 
into equation (-B4) and taking the first partial derivative 
with respect to (3, we find that for any positive 5, dyi/d(3 
and 8112/ dp are always negative. Therefore, both yi and y2 
solution branches decrease with increasing /J. The |m| > 3 
cases can be proven in a similar manner. 

For spiral cases, wc follow the same procedure to find 
that for \m\ > 2 cases both j/i and j/2 solution branches 
decrease with increasing /? while for the special case of |m| = 
1, the j/i and y2 solution branches increase and decrease with 
increasing /J, respectively. 



APPENDIX C: 

One can readily show that / tjf decreases with increasing 
D3 for both aligned and unaligned cases when |m| > 2 and 

\m\ > 1, respectively. For aligned cases, we simply rearrange 
equation (53) into the following form of 

(m2-2)(l + 5) , 2(m2-l)(l + <5) 



(CI) 



|m| \m\{D1 + l) 

which shows obviously that fJp / decreases with increasing 
Dl when \m\ > 2. 

For spiral cases, we can also arrange equation (92) into 
the parallel form of 

(m2-2)(l + 5) 



Mn(a)(m2 + a2 + l/4) 



■a^ + l/l)](i^<5) 



(C2) 



A/'^(Q)(m2 + a2 + l/4)(Di + 1) 

which again shows clearly that yfi j ii" decreases with increas- 
ing _Ds when \m\ > 1 for any radial wavenumber a. Note in 
particular that when |m| — 1, the coefficient of l/{Dg + 1) 
on the right-hand side of equation (C2) is positive. 




Figure Dl. The complete solution curves of marginal stability 
for spiral cases in a composite disk system when m = 0, 5 = 0.2 
and P = 1.5. The negative branches are shown here explicitly on 
purpose. 

system. For the cases of marginal stability for spiral con- 
figurations with ?n = 0, wc thought that there might be 
something interesting from the additional solution of a com- 
posite SID system in contrast to the single SID case (Shu et 
al. 2000). It turns out that this additional solution of is 
negative by numerical exploration. The numerical example 
of Figure Dl shows the typical structure of solution curves 
for stationary logarithmic spiral configurations when m — 0, 
S = 0.2 and P = 1.5. For mathematical completeness, the 
negative solution branches are shown explicitly. In Fig. 5 of 
the main text, the negative solutions of are discarded. 



APPENDIX D: 

In general, for two similar dynamical systems coupled to- 
gether, we expect two sets of oscillation modes (e.g., Van- 
dervoort 1991a, b). This should be true for the current prob- 
lem of a composite SID system as compared to a single SID 
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